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PREFACE 

This  is  the  forty-eighth  in  a  series  of  reports  growing  out  of  the  study  of 
radar  cross  sections  at  The  Radiation  Laboratory  of  The  University  of  Michigan. 
Titles  of  the  reports  already  published  or  presently  in  process  of  publication  are 
listed  on  the  following  pages. 

When  the  study  was  first  begun,  the  primary  aim  was  to  show  that  radar 
cross  sections  can  be  determined  theoretically,  the  results  being  in  good  agreement 
with  experiment.  It  is  believed  that  by  and  large  this  aim  has  been  achieved. 

In  continuing  this  study,  the  objective  is  to  determine  means  for  computing 
the  radar  cross  section  of  objects  in  a  variety  of  different  environments.  This  has 
led  to  an  extension  of  the  investigation  to  include  not  only  the  standard  boundary- 
value  problems,  but  also  such  topics  as  the  emission  and  propagation  of  electro¬ 
magnetic  and  acoustic  waves,  and  phenomena  connected  with  ionized  media. 

Associated  with  the  theoretical  work  is  an  experimental  program  which 
embraces  (a)  measurement  of  antennas  and  radar  scatterers  in  order  to  verify  data 
determined  theoretically;  (b)  investigation  of  antenna  behavior  and  cross  section 
problems  not  amenable  to  theoretical  solution;  (c)  problems  associated  with  the 
design  and  development  of  microwave  absorbers;  and  (d)  low  and  high  density 
ionization  phenomena. 


R.  E.  Hiatt 
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I 

INTRODUCTION 

1.1  Preliminary  Remarks 

This  is  the  second  in  a  series  of  reports  aimed  at  summarizing  the  available 
information  about  the  electromagnetic  scattering  properties  of  selected  bodies  of 
simple  shape.  The  first  body  considered  was  the  sphere  (Goodrich  et  al,  1961), 
which  from  the  bibliographer’s  standpoint  must  certainly  be  the  most  popular  of  all 
scattering  shapes .  The  present  work  is  devoted  to  the  cone  which,  while  less  abun¬ 
dantly  represented  in  the  literature,  still  presented  a  serious  selection  problem.  In 
the  case  of  the  sphere  it  was  immediately  obvious  that  a  large  amount  of  related 
reference  material  would  have  to  be  omitted  or  discussed  only  briefly  if  a  document 
of  manageable  size  were  to  result.  With  the  cone,  the  temptation  to  deal  with  the 
subject  exhaustively  is  stronger,  but  still  it  was  felt  that  such  an  attempt  would 
jeopardize  the  objectives  of  the  program  which  are  to  present  in  readily  usable  form 
available  theoretical  and  experimental  results  in  sufficient  detail  to  satisfy  most 
engineering  needs  and  provide  a  useful  guide  to  those  who  wish  to  delve  more  deeply 
into  the  subject.  Hence  a  certain  amount  of  arbitrariness  was  exercised  not  only  in 
deciding  what  material  was  to  be  excluded  but  also  in  the  relative  emphasis  given 
those  topics  which  survived  this  editorial  surgery. 

The  present  work  is  mainly  concerned  with  the  scattering  of  time  harmonic 
electromagnetic  waves  by  a  perfectly  conducting  cone,  either  finite  or  infinite  in 
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extent,  embedded  in  a  homogeneous,  isotropic  medium.  The  corresponding  scalar 
problem  involving  scattering  of  acoustic  waves  by  "hard"  and  "soft"  cones  is  also 
considered.  The  arbitrary  but  prevalent  division  into  scattering  and  radiation  pro¬ 
blems,  depending  on  whether  the  source  of  energy  is  off  or  on  the  cone  surface,  is 
observed.  Thus  in  restricting  attention  to  scattering  problems,  a  large  amount  of 
work  (e.g.  Felsen,  1957c;  Bailin  and  Silver,  1956;  Goodrich  et  al,  1959)  on  the  radi¬ 
ation  from  various  slot  configurations  or  dipole  distributions  on  conical  surfaces  is 
omitted. 

The  choice  of  the  cone  as  the  second  shape  to  be  treated  in  this  series  was 
dictated  by  two  reasons.  In  the  first  place  the  infinite  cone  is  one  of  the  few  shapes 
for  which  an  exact  theoretical  treatment  is  possible.  Since  the  approximations  which 
are  Introduced  in  order  to  make  tra  :able  the  scattering  problem  for  more  compli¬ 
cated  bodies  either  arise  from  or  are  tested  in  those  few  cases  where  exact  solutions 
exist,  familiarity  with  these  cases  seems  essential  for  an  understanding  of  the 
modern  methods  of  handling  scattering  problems.  The  second  reason  for  choosing 
the  cone  lies  in  the  fact  that  the  finite  cone  is  a  possible  re-entry  vehicle  configura- 
tlon,  so  that  its  scattering  properties  are  of  more  than  casual  interest. 

Unlike  the  sphere  for  which  exact  as  well  as  approximate  theoretical  results 
can  be  compered  with  experiment,  the  infinite  cone,  for  which  an  exact  theory  is 
also  available,  is  clearly  an  untenable  experimental  shape;  on  the  other  hand,  the 
finite  cone,  for  which  a  considerable  amount  of  experimental  data  is  available,  has 
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thus  far  resisted  exact  treatment.  At  one  time  it  was  felt  that  results  for  the  infinite 
cone  would  be  useful  in  predicting  the  finite  cone  behavior.  While  this  is  still  true 
in  a  sense  (tip  scattering),  it  is  now  realised  that  the  most  vital  scattering  charac¬ 
teristic  of  the  finite  cone  is  its  termination,  which  makes  it  fundamentally  different 
from  the  infinite  cone.  Nevertheless,  in  keeping  with  the  historical  development  of 
the  problem,  both  finite  and  infinite  cones  will  be  treated  here. 

1.2  Historical  Survey 

Since  the  circular  cone  is  a  level  surface  (9  -  9q)  in  spherical  coordinates  it 
was  a  natural  step  from  the  development  of  spherical  harmonics  of  nonintegral  order 


FIGURE  1 


3 


THE  UNIVERSITY  OF  MICHIGAN 

3648-2-T 


in  the  latter  part  of  the  nineteenth  century  to  the  investigation  of  potential  problems 
for  the  cone.  Although  Green  claimed  to  have  found  the  potential  near  the  vertex  of 
a  conducting  cone  in  an  electrostatic  field  as  long  ago  as  1828  (Green,  1828;  Mac¬ 
donald,  1900a)  this  work  was  never  published  in  detail  and  it  was  not  until  much  later 
that  the  problem  was  adequately  treated  in  the  literature  (Mehler,  1870;  Heine,  1878; 
Hobson,  1889;  Macdonald,  1900a).  Work  on  the  dynamic  problem  began  more  re¬ 
cently.  Macdonald,  in  his  aptly  designated  prize  essay  (Macdonald,  1902),  discussed 
the  problem  of  a  perfectly  conducting  cone  excited  by  what  is  now  known  as  an  axial 
(or  vertical)  electric  dipole  (dipole  moment  oriented  along  the  axis  of  symmetry). 
With  remarkable  insight,  Macdonald  formulated  and  solved  the  problem  as  a  series 
of  spherical  harmonics  whose  coefficients  he  determined  from  the  known  solution  of 
the  potential  problem.  Since  he  never  worked  directly  with  the  Helmholtz  equation 
but  with  a  related  partial  differential  equation  which  he  derived  from  Maxwell's  equa¬ 
tions  he  apparently  was  unaware  of  (or  at  least  unconcerned  with)  the  fact  that  his 
work  also  held  the  solution  of  the  scalar  Dlrichlet  problem:  a  soft  cone  with  an 
acoustic  point  source  on  the  axis  of  symmetry.  The  fact  that  electromagnetic  scat¬ 
tering  problems  for  cones  and  spheres  with  axial  dipole  sources  can  be  formulated 
in  terms  of  Dlrichlet  and  Neumann  scalar  scattering  problems  follows  from  the 
representation  of  an  electromagnetic  field  in  terms  of  radial  Hertz  vectors  which 
are  closely  related  to  Debye  potentials.  This  work,  for  the  most  part  more  recent 
than  Macdonald's,  is  discussed  in  more  detail  in  the  next  section  where  one  of  the 


resulting  scalar  problems  is  treated. 
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The  solution  of  the  scalar  Dlriohlet  problem  was  presented  by  Carslaw  (1910) 
in  essentially  the  same  form  as  Macdonald's  although  the  applicability  of  Macdonald's 
result  was  apparently  unknown  to  him.  A  few  years  later  Carslaw  (1914)  published 
the  solution  of  the  soalar  Neumann  problem  for  the  Infinite  cone.  In  this  later  paper 
Carslaw  indicates  awareness  of  Macdonald's  work,  calling  attention  to  some  weak 
points  in  the  derivation  (though  the  result  is  correct).  The  objection  concerned  Mac¬ 
donald's  argument  that  the  coefficients  in  the  series  solution  of  the  wave  equation 
were  the  same  as  those  for  the  corresponding  potential  problem,  i.  e.  were  indepen¬ 
dent  of  k.  Subsequently  Macdonald  (1915)  demonstrated  the  validity  of  his  argument 
and  this  was  acknowledged  by  Carslaw  (1916).  Nevertheless  Carslaw's  treatment  is 
more  elegant,  if  not  more  rigorous,  and  has  since  served  as  the  classic  work  on  the 
problem  to  the  extent  that  Macdonald's  work  is  apparently  overlooked  by  modern 
writers  on  the  subject. 

Following  Carslaw's  thorough  treatment  of  the  scalar  problem,  it  was  not 
until  the  resurgence  of  interest  in  classical  scattering  problems  occasioned  by  the 
development  of  radar  that  the  cone  received  further  study.  In  a  series  of  reports, 
Hansen  and  Schlff  (1948)  discussed  a  number  of  scattering  problems  including  that  of 
finding  the  field  scattered  by  a  perfectly  conducting  cone  when  a  plane  electromag¬ 
netic  wave  is  incident  along  the  axis  of  symmetry.  This  is  fundamentally  a  vector 
problem  since  the  vector  plane  wave  has  no  simple  representation  in  terms  of  mean¬ 
ingful  scalar  sources.  They  constructed  their  solution  in  the  form  of  infinite  series 
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of  Legendre  functions  of  non- integral  order  using  Hansen's  vector  wave  functions 
(Stratton,  1941;  Mentzer,  1955). 

Evaluating  the  coefficients,  however,  involved  a  rather  questionable  use  of 
asymptotic  forms  which  makes  the  procedure  subject  to  criticism  even  though  the 
correct  results  are  found.  An  additional  difficulty  is  present  when  the  source  is  at 
infinity  (plane  wave  incidence),  rather  than  at  some  finite  distance  from  the  vertex, 
since  the  far  field  expressions  are  divergent  sums.  This  is  also  true  in  the  scalar 
case  for  plane  wave  incidence  (Siegel  and  Alperin,  1952).  Special  summation  tech¬ 
niques  were  introduced  by  Schensted  (1953)  to  obtain  meaningful  results  for  the  back- 
scattering  cross  section  in  both  vector  and  scalar  cases  for  large  and  small  cone 
angles.  (See  also  Siegel  et  al  1953a,  1955b).  These  cases  were  treated  indepen¬ 
dently  by  Felsen  (1953,  1955,  1957a,  b).  He  found  the  Green's  functions  for  scalar 
and  vector  problems  and  based  his  calculations  of  back  scattering  cross  section  for 
small  and  large  cone  angles  on  these  more  general  results.  The  calculations  agreed 
with  the  results  of  Schensted  (1953)  and  Siegel  et  al  (1953a,  1955b)  but  by  working 
with  integral  representations  of  the  radial  functions  which  occur,  Felsen  avoided  the 
necessity  of  handling  divergent  sums.  However,  when  the  quantities  of  interest  are 
the  field  components  rather  than  the  cross  section,  both  series  and  integral  expres¬ 
sions  have  drawbacks.  In  particular,  any  attempt  to  calculate  the  currents  on  the 
cone's  surface  meets  with  additional  convergence  problems  in  the  series  expression 
and  a  considerable  complication  of  the  integral  expression. 
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A  slightly  different  method  of  constructing  the  vector  Green's  function,  using 
the  Lorentz  reciprocity  theorem,  was  given  by  B&illn  and  Silver  (1956),  but  calcula¬ 
tions  were  carried  out  only  for  the  radiation  problem  when  the  source  Is  on  the  cone 
surface. 

The  above  results  are  given  explicitly  in  Section  3.  Although  approximations 
may  be  employed  in  evaluating  the  expressions  for  the  field  or  the  cross  section,  the 
formulae  are  exact  In  the  sense  that  they  are  obtained  from  the  unique  solution  of  a 
well-set  mathematical  problem,  and  any  particular  calculation  can  be  performed  as 
accurately  as  time,  money  and  patience  will  permit.  This  is  borne  out  in  the  recent 
work  of  Goryanov  (1961)  who  reports  on  extensive  calculations  of  the  scattered  far 
field  when  a  plane  electromagnetic  wave  is  Incident  along  the  axis  of  symmetry  of 
the  cone. 

The  preceding  paragraphs  just  about  complete  the  survey  of  existing  exact 
treatments  of  scattering  by  a  semi-infinite  circular  cone.  A  more  general  problem, 
scattering  by  an  elliptic  cone  (of  which  the  circular  cone  is  a  special  case),  has  been 
solved  by  Kraus  and  Levine  (1961).  Although  of  special  interest  in  the  limiting  case 
where  the  cone  is  a  plane  angular  sector,  their  results  are  not  particularly  appro¬ 
priate  for  the  circular  cone  due  to  the  extreme  complexity  of  the  series  of  products 
of  Lam£  functions  in  terms  of  which  their  solution  is  expressed. 

Exact  treatments  of  the  finite  cone  problem  are  few  indeed.  Although  there 
are  uncountably  many  ways  in  whioh  the  cone  can  be  terminated,  what  little  exact 
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work  that  does  exist  deals  almost  exclusively  with  the  spherically  capped  right  cir¬ 
cular  cone+.  This  configuration  is  the  intersection  of  two  coordinate  surfaces  of  the 
same  spherical  coordinate  system  (r  =  a  and  6  -  Bq  in  Fig.  2).  Since  the  wave  equa¬ 
tion  is  separable  and  solutions  exist  for  each  of  the  surfaces  individually,  it  is  pos¬ 
sible  to  divide  space  into  two  distinct  regions  (see  Fig.  2)  in  which  different  series 


/  \ 

/  \ 


FIGURE  2 


representations  of  the  field  are  employed  depending  on  the  conditions  at  the  bound¬ 
aries  of  the  regions.  Thus  in  I  the  series  is  chosen  to  satisfy  a  condition  on  the 

^Referred  to  in  some  circles  as  the  "spherical  sector". 
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cone  whereas  in  Q  the  representation  must  satisfy  a  condition  on  the  relevant  portion 
of  the  surface  of  the  sphere.  On  that  part  of  the  boundary  of  the  two  regions  which 
has  no  physical  significance,  the  two  representations  must  be  equivalent.  This  mani¬ 
fests  itself  as  a  condition  sufficient  to  determine  an  infinite  set  of  constants,  but  the 
explicit  determination  usually  requires  the  solution  of  an  Infinite  system  of  linear 
equations.  This  procedure  was  followed  by  Northover  (1962)  who  considered  the 
problem  of  the  scattering  of  electromagnetic  energy  when  an  electric  dipole  is  ori¬ 
ented  along  the  axis  of  symmetry  of  a  perfectly  conducting  spherically-capped  cone. 

A  similar  procedure  was  adopted  by  Rogers  et  al  (1962)  in  considering  the  problem 
of  a  plane  electromagnetic  wave  incident  along  the  axis  of  symmetry.  Since  the  plane 
wave  is  the  limit  of  a  transverse  dipole,  this  problem  is  essentially  different  from 
the  one  treated  by  Northover.  The  big  drawback  in  both  cases  is  the  necessity  of 
dealing  with  infinite  systems  of  equations.  A  possibly  fruitful  way  to  handle  this 
problem  was  proposed  by  Plonus  (1961,  1962)  in  his  treatment  of  the  closely  related 
problem  of  calculating  the  radiation  pattern  of  a  blconical  antenna. 

A  technique  for  finding  the  exact  solution  of  the  scalar  wave  equation  satis¬ 
fying  Dirlchlet  boundary  conditions  on  the  intersection  of  two  surfaces  in  terms  of 
the  solutions  of  the  corresponding  Dirlchlet  and  Neumann  potential  problems  for  each 
surface  has  been  developed  by  Darling  (1960)  and  applied  to  the  spherically-capped 
cone.  The  method  is  restricted  to  low  frequencies  since  the  solution  is  expressed 
in  terms  of  an  infinite  series  which  converges  only  for  frequencies  in  a  bounded 
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range  (0  <  k  <  kQ).  Within  the  radius  of  convergence,  however,  the  series  repre¬ 
sents  the  exact  solution. 

Unfortunately,  these  exact  treatments  of  the  finite  cone  problem  do  not  con¬ 
stitute  a  sufficient  base  from  which  significant  quantitative  results  can  be  obtained. 
Since  the  finite  cone,  with  various  base  terminations,  is  an  object  whose  scattering 
properties  are  of  practical  importance,  approximation  methods  have  been  developed 
and  experiments  performed  in  order  to  obtain  the  desired  numerical  results. 

Experiments  are  limited  in  that  they  yield  results  at  particular  frequencies 
for  particular  models.  The  sensitivity  of  the  field  to  small  changes  in  frequency  or 
model  size  is  difficult  to  determine  since  it  is  impractical  (to  say  the  least)  to  try  to 
provide  a  continuum  of  models  or  operating  frequencies. 

Theoretical  approximations  are  also  limited,  though  most  often  to  a  range  of 
frequencies  rather  than  a  particular  frequency.  The  frequency  spectrum  is  divided 
roughly  into  three  regions.  The  low  frequency  end  of  the  spectrum  (wavelength  large 
with  respect  to  any  linear  dimension  of  the  scatterer)  is  usually  referred  to  as  the 
Rayleigh  region.  Clearly  this  has  meaning  only  for  finite  bodies,  for  present  pur¬ 
poses,  finite  cones.  Following  the  work  of  Rayleigh  (1897)  the  far  zone  scattered 
field  is  found  to  be  proportional  (to  a  first  approximation)  to  the  volume  of  the  scat¬ 
terer  regardless  of  its  detailed  geometric  configuration.  Some  refinements  of  this 
theory  as  well  as  explicit  results  for  the  scattering  cross  section  of  a  finite  cone 
have  been  given  by  Siegel  (1959)  and  further  extended  tty  Siegel  (1962,  1963). 
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The  high  frequency  end  of  the  spectrum  (wavelength  short  with  respect  to  the 
linear  dimensions  and  radii  of  curvature  of  the  scalterer)  Is  generally  referred  to  as 
the  optics  region  since  it  is  In  this  limit  that  the  laws  of  geometric  optics  are  appli¬ 
cable  (see,  for  example.  Born  and  Wolf,  1959).  For  many  purposes,  geometric  op¬ 
tics  Is  too  gross  an  approximation  and  refinements  are  sought.  One  such  Is  the 
physical  optics  approximation,  which  can  be  obtained  from  the  integral  representa¬ 
tion  of  the  field  everywhere  In  space  In  terms  of  the  field  on  the  scatterer  by  using 
the  geometric  optics  approximation  of  the  field  on  the  surface  of  the  scatterer.  This 
was  applied  to  the  cone  by  Spencer  (1951)  who  attributes  some  of  the  work  to  P.  M. 
Austin.  Hansen  and  Schiff  (1948)  cite  an  earlier  version  of  Spencer's  work  which 
also  credits  Crout,  who  worked  with  Austin  at  the  M.  I.  T.  Radiation  Laboratory.  A 
general  discussion  of  geometric  and  physical  optics  as  applied  to  the  computation  of 
back  scattering  (monostatic)  radar  cross  sections  is  given  by  Kerr  (1951).  The 
physical  optics  result  for  back  scattering  from  a  semi-infinite  cone  viewed  noee-on 
Is  given  by  Siegel  and  Alperln  (1952).  Although  the  applicability  of  physical  optics 
when  the  scattering  object  has  Infinite  length  Is  subject  to  considerable  question,  the 
result  was  shown  to  be  remarkably  good  when  compared  both  with  more  exact  theor¬ 
etical  results  (Siegel  et  al,  1955b;  Felsen,  1955)  and  experimental  data  (Sletten, 
1952).  More  general  results  are  given  by  Siegel  et  al  (1955a)  where  bistatic  physical 
optics  cross  sections  are  given  for  both  finite  and  Infinite  cones  when  the  transmitter 
is  on  the  axis. 
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Id  the  case  of  finite  cones,  the  physical  optics  approximation  is  least  reliable 
in  determining  the  effect  of  the  base,  and  since  this  is  usually  the  dominant  factor  in 
determining  what  the  field  will  be,  the  need  for  an  Improvement  over  physical  optics 
is  clear.  This  is  provided  for  the  flat- backed  finite  cone  by  the  geometrical  theory  of 
diffraction  (or  modified  geometric  optics)  of  Keller  (1957,  1958,  1960)  which  takes 
into  account  not  only  reflected  rays  but  diffracted  rays  as  well.  Expressions  for  the 
back  scattered  field  and  cross  section  as  a  function  of  angle  of  incidence  (plane  wave 
excitation)  are  given  by  Keller  (1960)  for  both  Dirichlet  and  Neumann  boundary  con¬ 
ditions  on  the  cone  in  the  scalar  case,  and  for  the  perfectly  conducting  cone  in  the 
electromagnetic  case.  This  last  result,  in  the  case  of  nose-on  incidence,  was  ob¬ 
tained  independently  by  Siegel  (1959)  and  Siegel  et  al  (1959). 

Keller  also  gives  results  for  a  cone  with  a  smoothly  joined  spherical  base 
(the  "cone-sphere"),  but  these  are  not  substantiated  by  experiment.  In  large  measure 
the  disagreement  is  attributable  to  Keller's  omission  of  any  contribution  from  the 
junction  of  the  cone  and  Bphere,  where  the  radius  of  curvature  is  discontinuous.  It 
is  new  believed  that  physical  optics  gives  a  reasonable  estimate  of  the  direct  return 
from  this  region,  but  it  is  also  true  that  no  canonical  problem  which  really  bears  on 
this  question  has  yet  been  solved. 

Between  the  Bayleigh  and  optics  regions  lies  the  resonance  region,  where  the 
wavelength  is  comparable  to  some  linear  dimension  of  the  scatterer.  It  is  this  re¬ 
gion  that  offers  the  most  difficulty  to  approximation  methods.  One  approach  to  the 
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problem  of  finding  approximations  which  are  useful  here  is  to  refine  the  small  and 
large  wavelength  approximations  and  effectively  narrow  the  width  of  the  resonance 
region.  More  terms  can  be  calculated  in  both  the  quasi-static  series  (cf  Senior  and 
Darling,  1963)  and  the  modified  optics  results  of  Keller  (1960)  but  when,  if  ever,  the 
regions  of  validity  overlap  is  not  known.  Another  attempt  to  bridge  the  resonance 
region  gap  based  on  physical  reasoning  is  described  by  Siegel  (1959)  and  Crispin  et 
al  (1963)  in  the  case  of  nose-on  back  scattering  from  flat-backed  finite  cones. 

Though  subject  to  the  limitations  previously  mentioned,  the  importance  of 
experimental  results  should  not  be  overlooked.  Thus  the  work  of  Sletten  (1952),  Olte 
and  Silver  (1959),  Keys  and  Prlmich  (1959a,  c,  d)  and  Blore  and  Royer  (1962)  has 
served  not  only  to  increase  knowledge  of  scattering  by  particular  cones  but  has  also 
played  a  valuable  role  in  the  development  of  approximate  theoretical  techniques. 

It  should  be  pointed  out  that  some  of  the  work  referred  to  contains  errors, 
some  typographical  and  some  more  serious.  Attention  is  called  to  these  errors  in 
the  sections  that  follow  and  formulae  presented  explicitly  have  been  corrected 
(wherever  possible)  and  converted  to  the  uniform  notation  used  throughout  this  re¬ 
port. 

In  concluding  this  section  the  authors  wish  to  extend  appreciation  to  the  many 
members  of  the  Radiation  Laboratory  who  assisted  in  the  preparation  of  this  report 
and  especially  to  J.  W.  Crispin,  Jr. ,  formerly  of  the  Radiation  Laboratory,  who  pre¬ 
pared  an  earlier  version. 
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EXACT  SOLUTIONS 

This  section  is  devoted  to  the  exact  formulation  and  solution  of  the  problems 
of  scattering  of  electromagnetic  energy  by  a  cone.  The  relation  between  vector  and 
scalar  boundary  value  problems  is  examined. 

In  some  cases,  a  single  scalar  function  can  be  used  to  solve  two  meaningful 
physical  problems,  one  scalar  and  one  vector.  One  such  case  is  treated  in  detail. 

2. 1  Precise  Formulation 

The  problem  with  which  this  report  is  concerned  is  that  of  finding  the  electric 
and  magnetic  field  vectors*  E  and  H  (and  quantities  derived  from  them)  external  to 
a  perfectly  conducting  cone  in  the  presence  of  various  incident  or  primary  fields. 

The  cone  is  assumed  embedded  in  a  homogeneous,  isotropic  and  perfectly  dielectric 
medium  of  permeability  p  and  permittivity  e,  which  medium  may  be  taken  as  free 
space.  The  propagation  constant,  k,  is  simply  related  to  n  and  e  as  follows: 

k  =  u  \fjje  =  y  (2.1) 

-iut 

where  A  is  the  wavelength  and  w  the  frequency.  The  harmonic  time  factor  e 
is  suppressed  and  MKS  units  are  employed  throughout.  At  all  ordinary  points  in 
space  the  behavior  of  the  field  quantities  is  governed  by  Maxwell's  equations 

+The  following  vector  notation  is  used  throughout:  vectors  of  arbitrary  magnitude 
will  be  underlined,  e.g.  E;  unit  vectors  will  be  denoted  by  carets,  e.g.  ?;  scalar 
products  indicated  by  dots,  e.g.  ?•  E;  and  vector  products  by  wedges,  e.g.  V^E. 
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VA  E  -  iuu  H  =  0  , 

(2.2) 

V,  H  +  iucE  =  0  , 

A” 

(2.3) 

•  E  =  V  S  -  0  . 

(2.4) 

The  homogeneous  equations  (2. 2),  (2. 3)  and  <2. 4)  do  not  describe  the  field  at 
source  points.  The  sources  treated  here  will  be  either  dipoles  or  plane  waves,  both 
of  which  are  defined  below.  The  presence  of  the  perfectly  conducting  cone  is  taken 
into  account  by  requiring  that  on  the  cone  the  tangential  component  of  £  must  vanish. 
With  the  addition  of  a  radiation  condition,  which  is  necessary  to  ensure  uniqueness, 
it  is  possible  to  formulate  a  well-set  boundary  value  problem  directly  in  terms  of  the 
electric  field.  It  is  customary,  however,  to  introduce  auxiliary  functions  from  whict 
the  field  quantities  may  be  derived,  and  to  formulate  the  problem  in  terms  of  these 
new  functions.  A  very  natural  way  of  introducing  both  the  Hertz  vectors  and  Debye 
potentials  stems  from  the  derivation  of  the  field  quantities  E  and  H  in  terms  of  a 
vector  potential  A  (see  Stratton,  1941).  If  A  and  f  are  vector  and  scalar  functions 
of  position  such  that 

fi72  +  k2)A=Vf  (2.5) 

then  E  and  H  can  be  defined  in  terms  of  A  by 

H  =  -iweV.A  ,  (2.6) 

“  A" 

s--£  v/ih*vaV-  ^ 
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It  Is  simple  to  show  that  E  and  H  so  defined  satisfy  Maxwell's  equations.  The 
divergence  condition  (2.4)  is  satisfied  because  the  divergence  of  a  curl  is  identically 
zero.  Equation  (2.3)  is  clearly  satisfied  since  it  is  taken  as  the  defining  equation 
for  E,  (2.7).  The  remaining  equation.  (2.2),  states,  after  substitution  of  (2. 1), 
(2.6)  and  (2.7)  that 

VA*.°.  (2.8) 

This  can  be  rewritten,  making  use  of  well  known  vector  identities,  as 

V(V2  +  k2)A  =  0  (2.9) 

A 

which  is  satisfied  by  virtue  of  equation  (2.5). 

Any  electromagnetic  field  can  be  so  derived;  that  is,  there  exists  a  vector 
potential  (in  fact  an  equivalence  class  of  potentials)  for  any  field.  In  particular, 
there  is  a  vector  potential  from  which  the  field  exterior  to  a  conducting  cone  can  be 
obtained.  However,  rather  than  investigate  the  restrictions  imposed  on  this  poten¬ 
tial  function  by  regularity,  single  valuedness,  and  other  physically-based  conditions 
on  the  field  quantities,  it  is  more  customary  to  restrict  the  form  of  the  potential  and 
phrase  the  problem  in  terms  of  this  smaller  class  of  potentials.  Of  special  interest 
is  the  construction  of  vector  potentials  from  scalar  potentials.  In  this  regard  if  <p 
is  a  scalar  solution  of  the  wave  equation 

(V2  +  k2)^  =  0  (2.10) 
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and  the  vector  potential  A  is  restricted  to  be  of  the  form 

A  =  ipa  (2.11) 

where  a  is  a  vector  independent  of  V  (but  may  be  a  function  of  position),  then  as 
shown  by  Senior  (I960),  a  must  either  be  a  constant  or  radial  vector.  (This  seem¬ 
ingly  preferred  status  enjoyed  by  rectangular  and  spherical  coordinate  systems 
results  from  the  restrictions  placed  on  A  and  does  not  pose  any  relativistic  dilem¬ 
mas.) 

With  the  potential  constructed  in  this  way,  the  class  of  electromagnetic  fields 
that  can  be  obtained  is  limited,  but  by  introducing  a  second,  independent,  vector 
potential  subject  to  the  same  restrictions  as  the  first  it  is  possible  to  derive  any 
electromagnetic  field  from  these  two  potentials  as  follows: 


E  =  A,  + 


(2.12) 


where  either 


or 


a)  c 

2  2 


b)  Aj  =  r  . 
2  2 


In  case  a)  c  it  a  constant  vector  (and  may  be  taken  to  be  of  unit  magnitude  by 

redefining  J^),  and  j)'  are  scalar  wave  functions,  1.  e.  satisfying  (2.10),  and  the 
2 
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A  and  so  determined  are  electric  and  magnetic  Hertz  vectors,  usually  denoted 
by  jje  and  7Jm  respectively.  If  a  rectangular  cartesian  coordinate  system  is  ori¬ 
ented  so  that  c  lies  along  the  z  axis,  i.  e.  c  =  ?  ,  then  the  rectangular  field  com- 
ponents  defined  by  (2.12)  are 


E  * 

y 


E  * 
z 


d% 

8i2 

H  = 

X 

dxdz 

+1“‘ir 

dxdz 

•te  9y 

*2ii 

H  = 

y 

*1 

8y8z 

dydz 

♦“IT 

d\ 

9 

8z2 

+  k2#, 

H  - 
z 

8z2 

♦fll 

(2.13) 


In  case  b),  r  is  the  radius  vector  r  *  3d  +yf  +z?  ,  and  i).  and  <pn  are 
—  —  x  '  y  z  14 

again  scalar  wave  functions.  The  functions  ^  and  ^  are  called  Debye  potentials 
and  the  A  and  A_  so  defined  are  often  referred  to  as  radial  Hertz  vectors.  The 
spherical  polar  field  components  defined  in  (2.12)  are 


E  =  ~(r^)  +  k2r^ 
r  3r2  1  1 


H  =  -^(r^2)  +  k%2 

3r 


F  .  aiH>  ,  m.  *2 

9  =  r  drdO  Bin  9~W 


1  a2(rV 


it*  *1 


9  r  drdd  sin  9  df 


_ _ i  a2(r*'i) 

*  rsind  drdji 


!  8  (r^2)  8^ 


H0  *  rein 8  8r80  +  **  99 


(2.14) 
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Note  that  If  t//  =  0  then  E  =  0,  i.  e.  the  electric  field  is  transverse  to  any  radius 
1  r 

vector  (TE  case),  whereas  if  0  then  H  =0,  i.e.  the  magnetic  field  is  trans- 

2  r 

verse  (TM  case). 

Actually  the  statement  that  any  electromagnetic  field  can  be  written  in  terms 
of  two  scalar  functions  using  (2. 12)  with  a)  or  b)  is  a  bit  too  strong.  Wilcox  (1957) 
proves  the  validity  of  the  Debye  potential  representation  of  any  field  defined  every¬ 
where  in  a  region  between  two  concentric  spheres.  Bouwkamp  and  Casimir  (1954) 
establish  a  similar  result  for  regions  exterior  to  a  sphere  containing  all  currents 
and  demonstrate  the  equivalence  of  the  Hertz  vector  and  Debye  potential  represen¬ 
tations  as  well  (cf  Sommerfeld,  1935).  Nisbet  (1955)  shows  that  ordinary  points  with¬ 
in  the  sphere  containing  sources  admit  of  a  representation  by  two  scalar  functions 
but  not  necessarily  Debye  potentials. 

In  cases  which  are  not  included  in  these  rigorous  treatments,  it  is  usually 
assumed  that  such  a  representation  is  valid  and  the  justification  is  found  in  the  rea¬ 
sonableness  of  the  consequences.  Thus,  in  the  case  of  the  field  scattered  by  an  infi¬ 
nite  cone  (a  body  which  cannot  be  enclosed  in  a  sphere  however  large)  Debye  poten¬ 
tials  still  play  an  important  role. 

In  formulating  the  boundary  value  problem  in  terms  of  scalar  wave  functions, 
it  is  clear  from  (2. 14)  that  requiring  either 

^2 

01  =  °,  and  —  =0 

0  =  0 
o 
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or 


02=°. 


and 


e=e 


Is  sufficient  to  ensure  that  the  boundary  conditions  =  0  at  6-  9q  are  satis¬ 

fied.  Wilcox  (1957)  proves  that  a  condition  sufficient  to  guarantee  that  E  and  H 
satisfy  the  Silver- Mtiller  radiation  condition,  viz. 


OoE=OoH  =  0  uniformly  in  r 


where 


0  = 


lim 


is  that  E  and  H  be  formed  from  Debye  potentials  ^  and  ^  38  shown  in  (2. 14)  and 
that  these  satisfy  the  scalar  Sommerfeld  radiation  condition 


lim 


00  r  \8?  - 


l)  *° 


uniformly  in  r. 


In  the  usual  formulation  of  a  scattering  problem,  the  total  field  is  considered 
as  consisting  of  two  terms,  an  incident  or  primary  field  (i.  e.  the  field  which  would 
exist  were  there  no  scattering  object  present)  and  a  scattered  field,  which  can  be 
thought  of  as  a  perturbation  or  correction  term  accounting  for  the  presence  of  the 
scatterer.  It  Is  therefore  necessary  to  discuss  the  kinds  of  primary  fields  which  are 
of  interest.  In  most  radar  scattering  problems,  the  target  is  assumed  to  be  illum¬ 
inated  by  an  Incoming  plane  wave  with  E  and  H  transverse  to  the  direction  of  propa- 
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gatlon.  In  the  mathematical  formulation,  however,  it  is  often  preferable  to  treat 
sources  at  a  finite  distance  from  the  scatterer  even  though  it  involves  introducing 
singularities  (of  a  prescribed  nature)  into  the  field  expressions.  This  is  particu¬ 
larly  true  when  the  scatterer  is  infinite  in  extent,  as  is  the  circular  cone  treated  in 
the  next  section.  In  the  first  place,  the  radiation  condition  is  imposed  uniformly  on 
all  field  quantities  when  all  sources  are  a  finite  distance  from  the  origin,  whereas 
with  plane  wave  excitation  only  the  scattered  field  or  perturbation  term  obeys  the 
radiation  condition,  and  if  the  scatterer  is  not  finite  in  extent,  this  term  must  be  fur¬ 
ther  separated  into  a  reflected  term  which  does  not  satisfy  the  radiation  condition 
and  a  diffracted  term  which  does.  Secondly,  by  considering  sources  at  finite  dis¬ 
tances,  it  is  sometimes  possible  with  the  help  of  Debye  potentials,  to  solve  a  mean¬ 
ingful  scalar  (acoustic)  scattering  problem  and  a  vector  (electromagnetic)  scattering 
problem  at  the  same  time.  This  is  not  true  with  plane  wave  incidence  and  the 
attractive  prosp  -  of  solving  two  problems  with  one  function  is  not  to  be  overlooked. 
Lastly,  with  the  general  solution  of  a  scattering  problem  for  an  arbitrary  dipole 
source  (essentially  finding  the  Green's  functions  for  the  particular  scatterer)  it  is 
possible  by  superposition  to  construct  solutions  for  any  source  distribution  and,  by 
passing  to  a  limit,  the  plane  wave  case  can  also  be  derived.  Of  the  various  sources 
only  dipoles  will  be  considered  since  more  complicated  sources  can  be  represented 
as  a  distribution  of  dipoles. 
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Let  (r,  6, 0)  and  (r^  (K ,  f-.)  denote  the  field  and  source  points  respectively 
and  R  the  distance  between  them  (see  Fig.  3).  A  scalar  source  at  the  point 
1*1 ,  By  is  then  given  by 

ikR 

**V  ■ 

We  note  in  passing  that  tp  can  also  represent  an  acoustic  velocity  potential. 


z 


FIGURE  3 
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An  electromagnetic  source  at  the  point  r^,  0^,  0^  can  be  derived  from  ip  in  a 
variety  of  ways: 

ikR 

a)  In  equation  (2. 13),  let  and  JL  -  0,  or  alternatively,  JT.  =  0, 

ikR  1  R  2  1 

_  g 

JT  =  .  In  the  first  case  the  field  components  given  by  (2.13)  are  those  of  a 

2  R  eikR  a 

Hertzian  electric  dipole  at  (r, ,  0, ,  0, )  with  TT  =  1  i  =  €.  i  .  In  the  second  case 

111  *“6  R  z  •‘lz 

the  equations  (2. 13)  define  the  field  components  of  a  Hertzian  magnetic  dipole  at 

ikR 

(r^,  0^)  with  T  Expressed  in  the  notation  of  Stratton  (1941) 

Ue  and  TJ^  are  electric  and  magnetic  dipoles  of  moment  ~  4jre?z  and 

m^=  ,  respectively, 

z 

b)  In  similar  fashion,  the  equations  (2. 14)  can  be  used  to  define  electromag- 

IkR  ikR 

6  6 

netic  sources  by  taking  either  ip  -  — and  ip  =  0  or  ip  =  0  and  ip  =  ■■  --  . 

1  «  6  1  2  It 

The  equations  (2. 14)  then  represent  the  field  components  of  a  radial  Hertzian  elec¬ 
tric  dipole  (first  case)  or  a  radial  Hertzian  magnetic  dipole  (second  case).  At  first 
sight  it  would  appear  that  in  order  to  treat  these  dipoles  in  the  same  manner  as  the 
dipoles  described  in  a)  (which  will  be  designated  axial  to  distinguish  them)  it  is 
necessary  to  define  the  dipole  moments  as 

2^  -  4xer  and  m^  =  4xr  .  (2.15) 

One  would  be  hard  put  to  give  a  physical  interpretation  to  a  dipole  moment 
which  varied  with  field  point.  However  it  is  a  relatively  simple  matter  to  show  that 
if  the  constant  vector  r^  replaces  the  variable  vector  r  in  the  definitions  of  dipole 


23 


THE  UNIVERSITY  OF  MICHIGAN 

3648-2-T 


moment  there  is  no  change  in  the  field  quantities.  That  is. 


/  ikR\  /  ikRN 


(2.16) 


and,  of  course,  this  remains  true  if  one  operates  on  both  sides  with  a  curl.  To 
verify  (2. 16),  we  first  note  that 


B  -  U-x,| 


which  implies  VR  = 


i-^il 


( V  operates  on  field  variables) 

and  also  ^  r  =  r,  =0.  With  these  facts  in  mind  it  is  clear  that 

ikR  ikR  .  /  ikRv 


va~  v1 


IKK  IKK  ,  /  1KK\ 

^(r- r2)  R  *  A  ^'-l*  =  dR  (  R  )  VRA  -  "~1 

dR  \  R  )  jr-  rj  A  ^  '  "1  °  * 


Hence  (2. 16)  is  an  identity  and  the  dipole  moments  can  be  defined  as 


(1)  a  .  U)  „ 

j>  =  4ver1  and  m  =  4rr^ 

which  are  truly  radial  dipoles  at  (r^ ,  0^  ^)  though  admittedly  normalized  in  a  pecu¬ 
liar  way. 

Because  of  the  convenience  of  the  Debye  potentials  or  radial  Hertzian  vectors, 
scattering  problems  with  radial  Hertzian  dipole  sources  are  often  encountered  in  the 
literature. 
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In  one  special  case,  of  course,  the  relation  between  radial  and  axial  dipoles 
is  particularly  simple.  This  occurs  when  the  source  is  on  the  z-axis  (0^  =  0)  but  not 
at  the  origin.  Since  the  rectangular  coordinate  syctem  employed  here  is  oriented  so 
that  the  z-axis  is  co-linear  with  the  axial  Hertz  vectors,  we  then  have  r,  =  r,  i  . 
Consequently,  the  field  components  defined  by  (2. 13)  when 


ik  Vx2+  y2+  (z -fj)2 

TJe  =  tt  =  SJT’i  2  *z 

fx  +y  +(z-r1> 


(2.17) 


Hn'-V.-0 


are  the  same  as  those  defined  by  (2.14)  when 


*1  rjR 


*2  =  ° 


(2.18) 


The  factor  l/r  normalizes  the  dipole  moment  to  be  4xer  ,  or  in  this  case  4xet  . 

1  1  z 

(An  analogous  relation  holds  between  axial  and  radial  Hertz  magnetic  dipoles.) 


Explicitly 


H  ■  -iueye 
x 


H  =  fcjcxe 

y 


R  /_ik  _1_ 
\R2  "  R3 

,/jk  J_' 

\r2'r3/ 


„  ,  .  ikR  k  ^  31k  3 


„  ,  .  ikR  hi  .  3ik  3 

Vy,z-r‘)e  V*1T?, 


( Ik  k2[*2+y2]  -1 

Vr2  R3 


H  =  0  E  =  e 
z  z 


3ik(z-r  )2  3(z-rt)2 

- 4 - + - 5 

R’  R  , 
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>i  2  2  ,  .2 

where  R  =  i  x  +  y  +  (z- r^)  . 


In  terms  of  spherical  coordinates  these  become 


H  =  0 
r 


E  =  e^fr,  f  -  cos  0g) 
r  l 


Hfl  ~  0 


=  iue  r  sin  9  e 


(?-?) 


Eg  =  e  (-rein  0f  + sin  0g) 


Er 0 


(2.20) 


where 

.  k2r  3ikr  3r 

R3  R4  ‘  R5 

2 

k  rr, cos  0-2  3ikrr,  cos  0  3rr,  cos  6 

2ik  .  1  .1  1 

g  *  “a  - 3 -  - 4 -  "  - 5 -  • 

RZ  R3  R4  R5 

When  r^  -  0,  the  case  invariably  treated  in  standard  texts  (e.g.  Stratton,  1941),  it 
is  easily  seen  that  r  =  R  and  the  above  components  then  reduce  to 


H  =  0 
r 


E  =  -2  cos  9  e 
r 


'  r  r  ' 


Hg=  0 


=  iue  sin  6  e 


(?-?) 


Eg  =  -sinfl  e 


Er°- 


(2.21) 


There  is  some  difficulty  with  the  radial  dipole  formulation  when  r^  *  0  due  to 


the  presence  of  the  factor  r,  in  the  definition  of  the  dipole  moment.  The  radial  elec 
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trie  dipole  of  moment  4xer,  *  4xer,t  (9,  is  still  zero)  has  zero  moment  when 

“l  1  z  l 

r,  =  0  and  therefore  cannot  radiate,  i.  e.  the  field  components  are  all  zero.  If,  on 

the  other  hand,  the  dipole  moment  is  normalized  to  be  4xcf  ,  then  (see  equation 

z 

2. 18)  is  apparently  singular  when  =  0  for  all  r.  Of  course,  if  the  field  com¬ 
ponents  are  computed  first,  the  limit  r^  — »  0  presents  no  difficulty,  but  in  spite  of 
this  it  is  more  convenient  to  avoid  the  radial  dipole  formulation  for  sources  at  the 
origin.  In  many  important  examples,  indeed  all  the  problems  discussed  in  this  re¬ 
port,  the  source  will  not  lie  at  the  origin,  it  being  more  convenient  to  choose  the 
coordinate  system  so  that  the  scattering  surface  rather  than  the  source  is  simply 
described,  and  the  radial  dipole  formulation  no  longer  presents  a  difficulty. 

It  is  now  possible  to  formulate  simultaneously  in  terms  of  Debye  potentials 
both  vector  and  scalar  scattering  problems  involving  a  cone  in  the  presence  of  ele¬ 
mentary  excitations.  First  define  the  ranges  of  field  and  source  variables  to  be 


o 

A 

H 

u 

0  4  9,  <9 

1  0 

0  4^^2w 

r  >  0 

04  94  9 

0 

04$42r  . 

g 

If  a  single  valued  function  ^  of  all  six  variables,  twlce-differentiable  in 
each  except  (possibly)  at  the  boundary  6  =  9q,  can  be  found  such  that 

+  «  0  (2.22) 
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r^oo  r  ~  =  0  uniformly  in  r,  (2.23) 

and 

=  0  ,  (2.24) 

o 

then 

a)  the  velocity  potential  of  the  field  scattered  by  a  perfectly  soft  cone  in  the 

ikR 

q  8 

presence  of  a  source  _  1  is  given  by  <p  ; 

XV  1 

b)  the  electromagnetic  field  scattered  by  a  perfectly  conducting  cone  in  the 

presence  of  a  radial  electric  dipole  (with  moment  4irer  )  is  given  by  equations  (2. 14) 
ikR  1 

with  i/jx  =  0®  +  ,^2  =  0; 

and  c)  the  field  components  of  b)  in  the  particular  case  when  0^  =  0  also  repre¬ 
sent  the  electromagnetic  field  scattered  by  a  perfectly  conducting  cone  in  the  pre¬ 
sence  of  an  axial  electric  dipole  of  moment 

B 

If  a  function  ^  can  b®  found  which  fulfills  all  of  the  conditions  required  of 

g 

except  that  the  boundary  condition  (2 . 24)  is  replaced  by 

*  0  (2.25) 

0*0 
o 

then  d)  the  velocity  potential  of  the  field  scattered  by  a  perfectly  rigid  cone  in  the 

ikR 

A  fl 

presence  of  a  source  — —  is  given  by  \p  ; 

K  2 
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e)  the  electromagnetic  field  scattered  by  a  perfectly  conducting  cone  in  the 

presence  of  a  radial  magnetic  dipole  (with  moment  4jrr.)  is  given  by  equations  (2. 14) 

ikR 

with  ^  -  o,  i>2  -  ifrj  +  V"  ' 

and  f)  the  field  components  of  e)  in  the  particular  case  when  0^  =  0  also  repre¬ 
sent  the  electromagnetic  field  scattered  by  a  perfectly  conducting  cone  in  the  pre¬ 
sence  of  an  axial  magnetic  dipole  of  moment  4xr  1  . 

1  z 

B  8 

A  detailed  derivation  of  the  functions  ^  and  ^  is  carried  out  in  the  next 
section  for  the  cases  when  the  source  is  on  the  axis  of  symmetry  of  the  cone,  so  that 


This  section  is  brought  to  a  close  with  a  brief  word  on  the  important  case  of 
plane  wave  incidence.  The  plane  wave,  either  scalar  or  vector  can  be  thought  of  as 
a  limiting  case  of  a  point  source:  namely  as  a  point  source  removed  to  infinity  with 

all  field  quantities  appropriately  renormalized.  For  example,  consider  the  scalar 

eik*  i-g - 2 - 

source  ~r—  when  9.  =  0.  We  then  have  R  =  Vr  +  r,  -  2rr.  cos  9  ,  and  for  suffi- 
xv  1  XI 

ciently  large  values  of  r^ 


ikR 
e _ 

R 


/v 


ikr^  -  ikrcos  9 


e 


ikr, 

8  1 

Deleting  the  factor - ,  i.  e.  renormalizing  the  source  strength,  now  leaves  only 

-ikz  Fl 

e  ,  which  is  a  plane  wave  propagating  in  the  direction  of  the  negative  z-axls. 

M  _  ilCZ 

Unfortunately,  however,  letting  ^  -  e  and  ^  =  0  (or  ^  =  0,  ^  *  e  )  in 
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the  equations  (2. 14)  will  never  produce  a  plane  electromagnetic  wave,  and  conse¬ 
quently  for  plane  wave  incidence  the  scalar  and  vector  problems  must  be  treated 
separately.  Although  it  is  possible  to  find  two  scalar  functions  which,  through  the 
equations  (2. 14),  produce  a  vector  plane  wave,  these  functions  have  no  meaningful 
interpretation  in  the  scalar  sense,  since  to  produce  the  plane  wave  E  =  t^e-^2, 

H  *  e_ik2  “Bing  (2. 14),  it  is  necessary  that 


-ikz 

,  _  xe 
^1  , 2,  2  2. 


and 


*2  = 


-ikz 


.  /2,  2. 
ku*i(x  +  y  ) 


k  (x  +  y  ) 

which  exhibit  a  delta  function  behavior  at  all  points  on  the  z-axls. 

The  vector  plane  wave  is  the  limiting  case  of  a  dipole  when  the  dipole  goes  to 

infinity  in  a  direction  normal  to  its  moment.  Thus  the  above  example  of  a  plane  wave 

IkR 

propagating  down  the  z-axis  is  the  limit  of  an  electric  dipole  TT  *  ~~  t  where 

k2R  X 

the  limit  and  renormalization  are  as  in  the  scalar  case.  For  this  source,  and  indeed 
for  a  transverse  dipole  in  general,  there  is  no  meaningful  scalar  problem  which  can 
be  formulated  simultaneously.  Though  Debye  potentials  may  still  be  employed,  they 
have  meaning  only  together  with  the  equations  (2.14)  which  define  the  field  quantities, 
whereas  for  radial  dipoles  the  appropriate  Debye  potentials  are  meaningful  in  them¬ 
selves. 

The  problem  of  finding  the  electric  field  scattered  by  a  perfectly  conducting 
leone  in  the  presence  of  a  dipole  of  arbitrary  location  and  orientation  is  one  of  deter- 


30 


THE  UNIVERSITY  OF  MICHIGAN 

3648-2-T 


mining  a  single-valued,  non-singular,  twlce-dlfferentiable  vector  function  ES  such 


that 

(1)  (V2  +k2)ES  *  0 

(2)  V  •  E®  -  0 

<3)  w.t/'S  £B)  * lkrE8}  ‘  0 

(4)  (E8  +  ESl  =  0  -  (e!  +  Ebl 

r  r  I  0  =  0  *  *10=0 


(2.26) 


where  either 


or 


with 


i  (s**  \ 

E  =  ^  VA  2  J  (electric  dipole  of  moment  4irej>) 
i  (e**  \ 

E  =  luyj  ^  I  — ml  (magnetic  dipole  of  moment  4xm) 


The  most  important  cases  are  those  where  •  jg  =  0  or  r^  •  m  =  0  (transverse 

dipoles).  When  r, .  p  =  0  or  r, .  m  =  0  there  is  the  alternate  formulation  in 
-in  -1A  ~ 

terms  of  Debye  potentials  (equations  (2.22)  through  (2.25)),  and  any  arbitrary  mo¬ 
ment  can  be  decomposed  into  a  transverse  component  and  a  radial  component. 

2.2  Particular  Solution 

Consider  a  point  source  of  radiation  located  on  the  axis  of  a  semi- infinite 
right  circular  cone  at  a  distance  r^  from  the  vertex.  In  terms  of  spherical  polar 
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coordinates  (r,  0, 0)  with  origin  at  the  vertex,  the  cone  is  defined  by  the  equation 

0=0  >»/  2 
o  ' 

so  that  i-9  is  the  half  cone  angle,  and  the  source  is  assumed  situated  in  the  ex- 
o 

terlor  region. 

Attention  will  be  confined  to  the  scalar  problem  (and  appropriate  vector  prob¬ 
lems  as  described  in  the  previous  section)  in  which  the  boundary  condition  at  the  sur¬ 
face  of  the  cone  is  either  ip  =  0  or  -  0  where  ip  is  the  total  field,  and  in  the 
neighborhood  of  the  source  ip  has  the  prescribed  behavior 

ikR 
e _ 

kR  ' 

V~2  2 

r  +  r j  -  2rr^  cos  0  .  The  Dirichlet  condition,  ip  =  0,  is  appropriate 
to  both  an  acoustically  "soft"  surface  and  an  electromagnetically  perfectly  conducting 
surface  in  the  presence  of  an  electric  dipole  (axial  or  radial)  whilst  the  Neumann  con¬ 
dition  (-^  =  0)  corresponds  to  both  a  "hard"  acoustic  surface  and  a  perfectly  con- 

C7C7 

ducting  surface  in  the  presence  of  a  magnetic  dipole  (axial  or  radial).  The  task 
which  then  confronts  us  is  the  determination  of  a  solution  ip  of  the  scalar  wave  equa¬ 
tion  which  is  free  of  singularities  in  the  region  0  <  0<  0q  except  at  the  source  but 
including  the  cone  vertex,  which  satisfies  the  boundary  condition  on  the  cone  surface, 
and  which  has  the  appearance  of  an  outgoing  wave  at  infinity  (see  equations  (2. 22) 
through  (2.25)). 
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The  form  of  the  solution  is  suggested  by  the  representation  of  the  source 
function  in  cylindrical  polar  coordinates.  From  Stratton  (1 941),  we  have 


“  i  ^2  (2n+l)  jQ(kr ^  hn(kr>)  P^) 


(2.27) 


where  n  -  cos  0,  j  (x)  and  h  (x)  are  the  spherical  Bessel  and  Hankel  functions  of 
n  n 

the  first  kind  defined  in  terms  of  the  corresponding  cylinder  functions  by  the  equation 


Zn(x)  =  V2x  Zn+V2(x)‘ 


(2.28) 


r<  stands  for  r  when  rCr^  but  for  r^  when  r  >r^.  and  conversely  for  r^.  Fol¬ 
lowing  Carslaw  (1914)  the  sum  is  now  expressed  as  an  integral  over  a  contour  C  in 
the  complex  v  plane: 


kR  2 


ikR  x  \  (2H-l)jj/(kr<)hi/(kr>)Pi/(-#i) 


sini'T 


(2.29) 


and  because  of  the  behavior  of  the  integrand  for  large  |  v  |,  the  contour  can  be  taken 
as  a  parabolic  path  running  from  v  =  (1  +  ix)oo  to  v  -  (1  -  ix)oo,  where  x  is  a  real 
positive  number,  and  intersecting  the  real  v  axis  between  v  =  -l/2  and  v  =  0. 

The  above  representation  of  the  source  function  suggests  that  for  the  total 
(incident  plus  scattered)  field  we  write 


*"2 


<2iH-l )  Ji/(kr<)  h^kr^  ijn) 
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where  ijn)  is  a  linear  combination  of  P^(m)  and  chosen  so  as  to  satisfy  the 

various  boundary  conditions.  Taking  first  the  problem  of  the  hard  cone  *  0  at 

Ov 

the  surface)  we  observe  that  if 


dM  I/O 
0 


(2.30) 


the  boundary  conditldn  at  0  =  0q  is  satisfied  automatically.  In  addition  the  solution 
is  finite  in  0  <  6  £  0q  except  at  the  source,  where  it  has  the  desired  infinity,  satis¬ 
fies  the  scalar  wave  equation,  and  has  the  appearance  of  an  outgoing  wave  as 
r-)a>.  Inasmuch  as  all  the  conditions  are  now  satisfied,  the  solution  for  the  Neu¬ 
mann  problem  is  therefore 


rP(-ii)~Pti)-PW7-P(-<i) 

*/  dM  vo  v  dM  v  o 

o  o 

►  dM 

J  Slni/T 

~  P  <M  ) 

dM  VO 

L.  0 

(2.31) 

This  solution  may  also  be  expressed  as  an  infinite  sum.  To  this  end  we  note 
that  since 

P  <-#x)  *  (-l)nP  (m) 
n  n 

the  bracketed  terms  in  the  integrand  ox  (2. 31)  are  zero  for  Integer  v,  and  conse¬ 
quently  the  only  singularities  of  the  integrand  are  poles  at  the  zeros  of  “■  P^q)  . 
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As  a  function  of  v,  these  are  real  and  distinct  (see  Carslaw,  1914)  and  are  non¬ 
integral  except  for  v  =  0.  This  last  is  a  zero  because  of  the  identity 

P  <4)s  1 
o 

and  since  v  =  0  is  also  a  zero  of  sini/r,  the  residue  at  this  point  requires  special 
consideration. 

From  the  definition  of  the  Legendre  function,  we  have 


P,<K)  =  ^(-v,  i*l; 


(see  for  example,  Magnus  and  Oberhettinger,  1949)  and  hence  for  small  v 


Vjv)  =  1  +  1/  log-^  +  CXi/2)  . 
The  residue  of  the  Integrand  in  (2. 31)  at  v  -  0  is  therefore 


i  1o(kr<lho(kr>1 '  TTT 


and  if  the  contour  is  now  closed  in  the  right  hand  half  plane,  we  have 


(2.32) 


(2.33) 


(2.34) 


~rr  P  .<-**> 


21  ,  „  „  „  ,  ,  -T- (2kfl  WW  *o  V 

-  *'  2_ - ^ - 3 - P»w 


■T~~  P  (M  ) 
dvdu  vo 
o 


(2.35) 


+We  remark  in  passing  that  Carslaw's  formula  corresponding  to  (2.35)  is  in  error 
by  virtue  of  his  failure  to  separate  out  the  residue  at  v  =  0,  though  his  result  cor¬ 
responding  to  (2.37)  is  correct. 
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where  the  summation  is  over  all  the  zeros  of  — —  P  (u  )  for  v>0. 

du  vo 
o 

The  above  expression  can  be  further  simplified  by  using  the  fact  that  for  v 


(i-n2) 


{p.w  £ 


P„M  '  V"’ 


ip«] 

dM  V  [  r 


sin  v* 


(Cars law,  1914).  If  v  is  a  zero  of  — —  P  to  ),  then 

du  i/o 
o 


(l-/i  ) P  to  ) T —  P  (-*<)=-—  sini/r 
o  i/  o  du  1/  o  x 
o 


which  can  be  substituted  into  equation  (2.35)  to  give 


I^-)o(kr<)ho<kr>)+21  (2^1)Ji/(kr<Jhi/(kr>) 


(i-m2)p  to  )-r~- p 

oi/o  di/du  1/  o 
o 


This  can  be  written  somewhat  more  compactly  as 


*2  “  2‘ 


(2.37 


(l-y2)P  to  )  — — P  to  ) 

u  V  Vo  di/dw  I/  O 

o 


where  the  summation  is  over  all  the  zeros  of  r —  P  (u  )  for  v>-\  2.  Since 

du  i/o 
o 

P  (u  )  =  p  ,(u  )  , 

1/  O  *  J  Q 


all  the  zeros  of  the  derivative  are  included  in  this  summation. 
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For  the  case  in  which  the  cone  is  soft  rather  than  hard,  the  treatment  is 
analogous  to  that  given  above.  Corresponding  to  (2. 31)  we  have 


(2i/+l)ji/(kr<)hj/(kr>) 

sini/ir 


-  y*  wi 

P  (u  )  [ 

I/O  J 


d v 


(2.38) 


which  can  be  written  as 


(2l'+1VkrA(kr>)  W  . . , 

\  ‘  z, - ^ - Trr.rvw 

v  ;W 


(2.39) 


(cf  equation  2.35),  where  the  summation  is  over  the  zeros  of  P^(uq)  for  i/  > -1/2. 
Since  i/  =  0  is  not  a  zero,  the  difficulty  which  was  found  for  the  hard  cone  does  not 
now  appear,  and  using  equation  (2. 36)  the  above  result  can  be  expressed  directly  as 


0, 


=  -2i  (2j/  +  l)ji/(kr<)hv(kr>) 


(l-p2)^-p  )^-P  ) 

o  dv  v  o  du  vo 
o 


(2.40) 


(cf  equation  2.37).  We  note  the  omission  of  a  factor  v  +  1/2  in  the  equations  cor¬ 
responding  to  (2.39)  and  (2.40)  in  Carslaw  (1910). 

With  all  the  above  formulae  computation  of  numerical  results  is  laborious  for 
two  reasons.  In  the  first  place  it  is  difficult  to  obtain  sufficiently  precise  values  for 
the  zeros  of  the  Legendre  functions  and  their  derivatives,  and  secondly  the  infinite 
series  converges  veiy  slowly  under  most  circumstances  of  practical  interest. 
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Deferring  until  later  a  discussion  of  the  first  point,  attention  is  now  focussed  on  the 
convergence  properties. 

That  the  series  actually  do  converge  can  be  seen  by  considering  the  behavior 
of  the  elements  of  the  summation  for  large  v.  Taking  first  the  radial  functions  we 


VW'2  1/7  3/2>J  {1+°<»2/‘'>}  <2-«> 

{' +  ,2  421 


valid  for  i/»x  .  Also 


P  \/ —  , — n  sin 

v  V  mt  Bin  0 


ji/+l/2)0+|j 


(see,  for  example,  Carslaw,  1914),  giving 


“V*1  ri — 
*  "tjii 


(2.43) 


O'+l^cos  ■j(l'+l/2)e  +  |l  (2.44) 


dvdw  y  .  3 
I  vi  sin  0 


(v +  1  /2)'0 sin  j(v+  1/2)0  +jJ,  .  (2.45) 


Hence,  in  equation  (2. 37)  the  terms  for  which  v  is  sufficiently  large  can  be  written 
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which  tends  to  zero  as  v  — »  oo,  and  does  so  with  a  rapidity  which  increases  as 

r '</r>  decreases .  A  glance  at  equation  (2 . 43)  reveals  that  the  zeros  of  P^(mq)  are 

approximated  by  r/Q  (m  -  l/4)  -  1/2  (m  integer)  wherever  (2.43)  obtains  (for 

0  dP  Ou  ) 

v,  hence  m,  large).  Similarly  the  zeros  of  — -  are  approximated  by 

<> 

f/Mm+  1/4)  -  l/2,  (m  integer)  when  v  (and  m)  are  large.  Substituting  this  last 


value  in  (2.46)  yields 


-i 

kr^ 


1 


9  sin©  (2m  +7:  -  ■^Q)  kr> 

O  i  T 


t  _1  mi 

.  .40  "2  0 

/kr<\  o  o 


sin  j(m+  1/4)  +  fj  . 


(2.47) 


mi/e 

/kr  \  o 


which  in  absolute  value  can  be  shown  to  be  less  than 


(- 

\Jm  \kr>J 


,  where  c 


is  independent  of  m  and  bounded,  so  that  by  the  comparison  test  the  series  (2.37) 

kr< 

appears  to  converge  absolutely  when  — —  <  1 .  A  similar  result  also  applies  to  the 

Kr> 

sums  in  equations  (2.35),  (2.39)  and  (2.40).  A  completely  rigorous  proof  of  conver¬ 
gence  would  also  show  that  this  convergence  is  unaltered  when  all  higher  order 
terms  in  l/v  in  the  expressions  for  the  Bessel  and  Legendre  functions  (2.41)  - 
(2.43)  are  included.  Such  a  proof,  to  the  best  of  the  authors'  knowledge,  has  never 
been  given. 

Assuming  that  these  higher  order  terms  are  truly  negligible,  a  study  of 
(2.46)  or  (2.47)  Indicates  that  the  convergence  is  improved  if  the  source  (or  the 
receiver)  is  moved  off  to  infinity,  and  since  this  case  has  considerable  practical 
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Interest  it  will  be  considered  next.  To  this  end,  let  us  assume  that  in  equation 
(2. 37)  the  sum  can  be  truncated  at  some  value  v  =  v^,  say.  For  sufficiently  large 
kr^  the  spherical  Hankel  function  can  then  be  replaced  by  its  asymptotic  expansion 
for  large  argument,  viz 


-i  1(kr>"l/2)  r 

v(kr>'  ~  6  (’ 


1  -  i 


vCv+D  .  \ 

kr>  ■■■} 


(2.48) 


and  if  only  the  first  term  is  employed,  the  expression  for  4j  becomes 

6 


ikr„ 

*2  =2'V 


-iv- 


5]  e  2  (2i/+l)  j^(kr^ 


P» 


(1-M2)P  (m  )-tt~P  (H  ) 
o  v  o  di/d M  V  o 
o 


(2.49) 


valid  for  kr>»  i/^i/^+l).  Providing  that  the  truncation  of  the  original  series  can  be 

Justified,  equation  (2.49)  gives  the  solution  for  either  a  receiver  at  a  large  distance 

from  the  vertex  or  the  point  source  at  a  large  distance.  For  the  limiting  case  of  the 

source  at  infinity  (plane  wave  incidence)  it  is  necessary  to  re-normalize  by  mul- 

ikr> 


tiplying  through  by  kr>e  before  taking  the  limit  kry-^oo.  This  is  equivalent 

-1  _ikr> 

to  deleting  the  factor  (kr^  e  in  (2.49). 


To  obtain  an  estimate  for  it  is  sufficient  to  apply  to  the  sum  in  (2.49)  the 

tame  analysis  used  in  regard  to  equation  (2.37).  For  v  sufficiently  large  that  the 
lessel  and  Legendre  functions  can  be  replaced  by  their  approximations  for  large 
order  (equations  2.41  through  2.45),  the  corresponding  term  in  (2.49)  then  becomes 
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W2  \l — ~  jkf<^2)  sin  {(v+l/2)6+*/j 
12 sin©  f'(v+  3/2)  6 


ir  e 


which  can  be  written  as 

v  ,  3 


e 


*"  2  +  2  CT  \V  .to  ((,.+  1/2)9 +  */<) 

Vsine  \2iH-3  /  2i/+3  9  W'  ' 


by  using  Stirling's  formula  for  the  gamma  function.  This  tends  to  zero  rapidly  as 
v  ->  co  and  to  judge  from  its  form  it  would  appear  that  a  sufficient  criterion  on  v 
is 


1 


V1»kr<- 


(2.51) 


Accordingly,  the  introduction  of  the  asymptotic  expansion  for  h  (kr^)  is  justified  if 

kr>»(kr<)2  ,  (2.52) 

and  in  the  limiting  case  of  plane  wave  excitation  the  total  field  is  (from  equation 
2.49) 


*r!I 


-iv|  Pv(m) 

e  1  (2i/+l)j^(kr)  - 2 - 

v  (i-m2)p  (m  )t~—  P  (M  ) 

o  v  o  di/dM  v  o 
0 


(2.53) 


where  v  is  such  that 


dP  (u  ) 

vo 


=  0.  The  corresponding  plane  wave  result  for  the 


soft  cone  (equation  2.40)  is 


«  -2 


E 


-ll/|  (2l/  +  l)j^(kr)P^((l) 

8  rww 

“-"o’  “ST"  IT- 


(2.54) 
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where  v  is  such  that  P  (u  )  =  0. 

v  o 


It  is  easily  demonstrated  (e.g.  Siegel  and  Alperin,  1952)  that  when  v  is  such 


dP 

that  (n  )  =  0,  then 
du  o 
o 


(1  -p2)  d2P  (m  ) 

_ 0_  V  o  .  . 

2i/+l  du  dv 


\  Km2*  <2-55> 


and  when  v  is  such  that  P  )  =  0.  then 

vo 


(l-/i  )  dP  <m  )  dP  (n  ) 

O  I/O  VO 


dM 


dv 


M2  *  =  tfkfi  \ 


2v  +  l 

°  ^1  ^1 

(2.56) 

If  these  expressions  are  used  in  (2. 53)  and  (2. 54)  the  total  field  for  both  boundary 
conditions  for  plane  wave  excitation  can  be  expressed  as 


ip  =  2  '  e 


-lv-  Jy(kr)PvW 


£  M2* 


(2.57) 


where  ip  -  ip,  (soft  cone)  with  the  summation  over  v  such  that  P  in  )  =  0  and 

dW  v  0 

ip  =  <p g  (hard  cone)  when  — — -  =  0.  Note  that  a  uniform  expression  can  also  be 

lobtained  for  the  point  source  solutions  (equations  2.37  and  2.40),  namely 


,  ■  2.  y  WW! 
*  ZM2* 


(2.58) 
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where  the  same  convention  Is  observed;  namely  <p  is  the  solution  of  the  Dirichlet 

problem  when  the  v'a  are  zeros  of  P^q)  and  0  the  solution  of  the  Neumann  pro- 

dP  <ti  ) 

blem  when  the  v's  are  the  zeros  of  — - — —  . 

*0 

The  solution  required  in  most  calculations  associated  with  the  cross  section 
is  given  by  (2. 57)  but  unfortunately  its  form  is  not  convenient.  In  the  first  place, 
the  sum  represents  the  total  field,  whereas  the  incident  and  scattered  fields  are 
individually  necessary  in  order  to  compute  a  cross  section,  and  the  separation  of 
(2.  57)  is  a  not  too  easy  task.  In  addition,  the  definition  of  the  cross  section,  viz 

lim  .  2  (s,,i  2 

o  -  4jt  r  di  Nj 

r  — *00  ' 

implies  that  the  receiver  is  placed  at  such  a  distance  from  the  body  that  the  scat¬ 
tered  field  has  acquired  the  characteristics  of  a  spherical  wavet  and  avoiding  the 
question  of  whether  this  definition  is  applicable  when  the  body  is  infinite  in  the  trans¬ 
verse  direction  (JBrysk,  1960)  the  consideration  of  equation  (2. 57)  for  large  values 
of  kr  is  fraught  with  difficulties. 

If  the  Bessel  function  j^(kr)  in  (2.53)  and  (2. 54)  is  replaced  by  the  first  term 
of  its  asymptotic  expansion  for  large  kr  the  resulting  series  diverge.  However,  in 
the  case  of  back  scattering,  Schensted  (1953)  and  Siegel  et  al  (1953a,  1955b)  use 
special  summation  methods  to  obtain  expressions  for  the  cross  section.  The  diffi¬ 
culty  inherent  in  treating  divergent  series  is  avoided  by  Felsen  (1955)  who  replaces 
the  Bessel  function  not  by  an  asymptotic  form  but  by  its  Sommerfeld  integral  repre- 
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sentation,  and  obtains  by  steepest  descents  analysis  an  approximate  evaluation  of 

both  the  field  and  the  cross  section  for  large  and  small  angle  cones.  The  results  of 

these  analyses  are  presented  in  Section  3. 

This  section  is  brought  to  a  close  with  a  brief  discussion  of  the  computational 

problem  of  finding  precise  values  for  the  zeros  of  the  Legendre  functions  and  their 

first  derivatives  (those  values  of  v  for  which  P  (cos  9)  or  —  P  (cos  9)  vanish  for 

v  da  v 

-  d  1 

fixed  values  of  9) .  In  addition,  the  zeros  of  — ■  P  (cos  9)  are  of  interest  in  the  vec- 

d0  v 

tor  problem  (see  section  3.1.3). 

For  large  values  of  v,  these  zeros  can  be  simply  approximated,  as  was  men¬ 
tioned  earlier  in  the  discussion  of  convergence.  However,  the  dominant  contribution 
to  the  series  (2.  57)  and  (2.  58)  (or  any  of  the  alternate  forms  given)  comes  from  the 
early  terms  and  consequently  the  small  zeros  (small  values  of  v  for  which 
P™(cos  9)  =  0)  are  quite  important  in  any  computation  of  the  field  or  the  cross  sec¬ 
tion. 

Obtaining  accurate  values  for  these  small  zeros  and  increasing  the  accuracy 

of  the  large  zero  approximation  are  not  trivial  problems.  Hobson  (1931)  summarizes 

the  work  of  Macdonald  (1900b)  and  Pal  (1918,  1919)  who  find  approximations  to  the 

zeros  as  a  function  of  argument.  These  expressions  were  employed  by  Carrus  and 

Treuenfels  (1951)  who  tabulated  the  first  50  zeros  of  P*(cos  9)  every  5°  for 

90°<  9<  175°  and  the  first  50  zeros  of  ~  P*(cos  9)  every  5°  for  90°<  9  <  130°. 

d0  v 
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Also  tabulated  were  the  Integrals 


\  0[p>>]2*‘ 


for  these  i/'s.  Errors  in  these 


tables  were  noted  by  Siegel  et  al  (1951).  For  the  special  value  6q=  165  the  first 
twenty  zeros  and  normalizing  integrals  were  recomputed  by  Siegel  et  al  (1952, 
1953b).  These  last  were  again  recomputed  with  even  more  accuracy  by  the  Institute 
of  Numerical  Analysis  (see  Siegel  et  al,  1953a). 

Even  if  the  zeros  were  tabulated,  however,  the  calculation  of  the  field  or 
cross  section  would  still  present  a  formidable  computational  problem.  Some  of  the 
difficulties  inherent  in  the  representation  as  infinite  sums  over  appropriate  zeros 
are  avoided  in  the  integral  representations.  Appropriate  integral  expressions,  if 
available,  should  be  used  in  any  extensive  computational  program  [jsee  Section  3 
and  Felsen  (1953,  1955,  1957a,  b,d,  1958)J  .  This  is  especially  true  in  the  light  of 
present  day  high  speed  computing  techniques,  as  is  borne  out  by  the  work  of  Gory- 
anov  (1961)  (see  section  3.  3). 
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in 

SUMMARY  OF  SEMI-INFINITE  CONE  FORMULAE 


This  section  presents  a  systematic  listing  of  existing  expressions  for  the  scat¬ 
tering  properties  of  semi-infinite  cones  with  either  plane  wave  or  point  source  exci¬ 
tation.  Various  approximations  are  given,  together  with  the  restrictions  on  the  para¬ 
meters  under  which  the  approximations  are  valid.  All  formulas  have  been  trans¬ 
formed  to  the  notation  employed  throughout  this  report.  Thus  ip 1  represents  a  scalar 

source  (either  point  source  or  plane  wave),  ip.  represents  the  field  scattered  by  a 

is  s 

soft  cone  (Dlrichlet  boundary  condition)  and  <p^=  <p  +< p^  is  the  total  field,  t p^  rep¬ 
resents  the  field  scattered  by  a  hard  cone  (Neumann  boundary  condition)  and 
1  s 

^  +^2  '  111  general  these  < p  will  be  functions  of  six  variables,  viz.  3  field  coor¬ 
dinates  and  3  source  coordinates.  When  ip  is  a  plane  wave,  the  number  is  reduced  to 
five  and  in  the  important  case  when  the  plane  wave  is  incident  along  the  axis  of  sym¬ 
metry,  the  field  is  a  function  only  of  the  three  field  coordinates. 

When  the  scattered  field  is  observed  far  from  the  vertex,  it  is  convenient  to 

ikr 


S  (0,0)  and  j  =  1  or  2.  For 


8  6 

introduce  the  scattering  amplitude  where  ip  ^  ^ 
general  plane  wave  incidence  the  scattering  amplitude  is  also  a  function  of  0  ,  0  . , 
the  angles  specifying  the  direction  of  incidence.  The  scalar  scattering  cross  sec¬ 
tions  will  be  denoted  by  a  where 
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In  general  a.  will  be  a  function  of  four  angular  variables,  but  most  of  the  results 
available  are  for  the  case  0^  =  0  (incidence  along  the  z-axis)  so  that  the  number  is 
reduced  to  two.  The  back  scattered  far  field  for  plane  wave  incidence  along  the  z- 
axis  is  an  oft -computed  quantity  and  will  be  denoted  by  Sj(0)  and  o^(0).  The  electro¬ 
magnetic  results  will  all  be  given  in  terms  of  the  electric  field  E ,  and  the  same 
superscript  convention  will  be  observed  as  described  above.  Thus,  E*  is  the  in- 

g 

cident  field,  E_  the  scattered  field,  and  S  the  far  zone  scattering  amplitude  (now  a 
vector),  a  (without  subscript)  will  denote  cross  section,  defined  as 


<t  = 


2 


andcr(O)  is  the  back  scattered  cross  section  for  incidence  along  the  z-axis. 

hi  the  presentation  of  formulae,  the  natural  division  into  vector  and  scalar 
problems,  with  the  latter  further  split  into  Dirlchlet  and  Neumann  boundary  con¬ 
ditions,  will  be  observed.  With  each  formula  will  appear  the  appropriate  conditions 
for  applicability. 
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3. 1  Dlrlchlet  Boundary  Condition. 

For  an  arbitrarily  located  scalar  source 


e _ 


in  the  presence  of  a 


soft  (4j\.  =  0)  semi -infinite  cone  of  half  angle  t  -9n  (equivalent  to  a  radial  elec- 

4r  er^ 

trie  dipole  of  moment  — - -  and  a  perfectly  conducting  cone)  the  total  field  is 


given  by  Felsen  (1957  a,  b)  as 

00  p 

emcoBm(0-01)  I  d/ 
m=o  J 


v 


(2v  +1) 


r  (v  +m+l  )}v  (kr^h^  (kr^ 
r  (v  -m+l)sin(v  -m)x 


P  m(cos  ©  ) 
V  o 


[p  m  (cos  0  )P  m  (-cos  ©)-P 
\_v  >  v  0  v 

%  V 


-m .  .  .  -m. 

(-cos 0  )P  (cos©,, 

>  V 


1) 


where 


=  Jv+]/2(x)>  Vx)=^722  h(^1/2(x), 

eo  =  lf  Sn=  2  for  m  =  1,  2,  3, . . . ,  + 

r  =  r  for  r  >  r^  , 


r<=  r  for  r  Cr^ 

=  r^  for  r^  <  r 


©  =  0  for  ©  <  9 ^ 

=  &  j  for  0  <9 


=  r^  for  r^  >  r 


9>  =  9  for  9  >  ©i 
=  ©  ^  for  ©  1  >  © 


and  the  contour  c  is  shown  in  Fig.  4. 


These  definitions  of  em  and  the  spherical  Bessel  functions  agree  with  those  in 
Magnus  and  Oberhettinger  (1949)  but  differ  from  those  in  Felsen  (1957b). 
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v  -plane 
FIGURE  4 


When  the  source  is  on  the  axis  of  symmetry  (0^  0)  equation  (3. 1)  simplifies 
considerably,  and  only  one  of  the  many  alternative  formB  for  the  resulting  egres¬ 
sion  is  given  here.  Alternate  expressions  can  be  found  in  Section  n  of  this  report, 
Felsen  (1955),  Carslaw  (1910)  and  Macdonald  (1902),  though  care  is  necessary  to 
make  sure  that  the  source  normalization  and  time  dependence  are  consistent  before 
comparisons  are  drawn.  An  additional  problem  is  present  in  the  older  references 
(Carslaw  and  Macdonald)  because  their  notation  is  usually  different  from  that  em¬ 
ployed  by  modern  writers;  e.g.  Macdonald's  (a,  0,7),  and  (x,  y,  z)  are  the  components 

~iv* 

of4*H  and  ~  E  respectively,  and  Carslaw' s  iMikr)  is  -jp  e  2  H^kr). 

ikR  g - ' 

The  total  field  for  a  source  6  ,  with  R  =  J  r  +  r^  -2  r  r^  cos  0  ,  is 
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^  =  2i 


j  (kr  )h  (kr  )P  (cos 0) 
v  <  v  >  v 


(3.2) 


where  HQ  =  cos  0Q  and  the  summation  is  over  the  positive  zeros  of  (mq),  that  is, 
over  the  set  of  v  such  that  v  >  -1/2  and  P^  (uQ)  =  0. 

A  symmetric  distribution  of  electric  dipoles  about  the  z-axis  is  discussed  by 
Macdonald  (1902)  and  a  simply  varying  circular  source  is  treated  by  Felsen  (1959). 

When  a  plane  wave  is  incident  from  an  arbitrary  direction  the  total  far  field 
is  given  by  Felsen  (1957b)  as 


-ikr|sin0  sin0icos(0  )+cos0 cosojj 


i  e 


ikr 


2kr{cos0  +COS0J  )log  — — 


(3.3) 


subject  to  the  restriction  kr  >>1,  0+0^<20q-t,  and  0^  r  (thin  cone). 

i  -ikz 

When  a  plane  wave  is  incident  along  the  axis  of  symmetry,  tp  =  e  ,  the 
total  field  is  given  in  Section  n  £also  Felsen  (1955)]  as 


-iv  r/2  j  (kr)P  (cos  0 ) 
e  '  *v  v _ 

•Cm2-" 


(3.4) 
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ilcz 

The  coefficients  of  -e  in  equations  (3. 8)  and  (3.  9)  are  plotted  in  Fig.  5  where  it  is 
seen  that  the  two  coincide  for  §^4  +;  hence  (3. 8)  must  be  used  for  £  <  4  while  the 
simpler  expression  (3.  9)  suffices  for  £  >  4  (Felsen,  1958). 

Corresponding  to  (3.  5)  and  (3.  7)  the  following  expressions  for  the  back  scat¬ 
tering  cross  section  obtain: 

x2 

0.(0)= - - - r-jr  (3.10) 

16*  (log  — ) 

*-B0 

when  0qa»  r,  and 


o^(0)= 


ie*(e0-f)4 


(3.11) 


when  0Q~  -  . 

The  cross  sections  given  in  (3. 10)  and  (3. 11)  are  plotted  in  Fig.  6  as  a  function  of  0Q. 


3.2  The  Neumann  Boundary  Condition. 


*1^1 


For  an  arbitrarily  located  scalar  source 


in  the  presence  of  a 


hard 


=  o) 


k\^i\ 

semi-infinite  cone  of  half-angle  #-0o  (equivalent  to  a  radial 

Q  ^  4- 

magnetic  dipole  of  moment  ~r  £  and  a  perfectly  conducting  cone)  the  total  field  is 

K  i- 

given  by  Carslaw  (1914)  £see  also  Felsen,  1957b]  as 


oo 

^2=2i  Y1  V081"^'^)^ 


(2v  +1)^  (kr. )hv  (kr>)P^  “(cos©  )P^  "Ws  0j) 


m=o 


(3.12) 


''The  amplitude  curves  coincide  at  £  **  2  but  the  phases  differ. 
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where  the  notation  is  the  same  as  defined  under  equation  (3. 1)  except  that  the  sum- 

dp-m 

mation  is  over  the  non-negative  zeros  of  •— —  (cos  6  ).  That  is,  the  set  of  v's  for 

d6°  cfr-m 

each  m  consists  of  those  v  for  which  v  >  -  1/2  and  -j—  (cos  6  )=0.  Note  that 

d9o  ° 

when  m=0,  v  =0  is  an  element  of  the  set. 

When  the  source  is  on  the  axis  of  symmetry  (0^  =0)  the  field  expression  is 
given  more  simply  (see  Section  II)  as 


,21  y  8 ) 

2  c—j  ><1 
v 


1 


(3.13) 


dp 


dP 

where  the  v  are  now  such  that  v  >  -1/2,  — — (cosfl  )=0.  Alternate  expressions  can 

d»0 

be  found  in  Carslaw  (1914)  and  Felsen  (1955,  1957a,  b). 

When  a  plane  wave  is  incident  from  an  arbitrary  direction  so  that 

j  -ikr(sin0  slnd^cos^  -^)+cos0cos0i] 

=  e 

the  scattered  field  is  given  by  Felsen  (1957b)  for  thin  cones  as 


,s  2 


(tt)  i/S;e  * 


(3.14) 
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where 


xf 


dx|x  - 5 — H'  vkr)(l/4+x^)p  ,(cos0)P  /(cose!) 

1  cosh'  xx  **  ix'V2 


/*°°  xf 

1=  dxJ  x  •-  --  ~~a — H^\kr)P*  ,  (cos0)P* 

1  J  cosh2,x  lx  ix-^  ix‘ 

-00 


,  (cos  01 ) 
Va 


0^»,  0+0,  <20„-»  . 

o  1  o 


In  the  far  zon^  equation  (3. 14)  can  be  slmplifed  (Felsen  1957b)  to  read 


i(r-0  )  ik: 

,s  o  e 

- 3  — 

2(cos  0  +  cos  0^ )  kr 


— —  M.+COS0  008  0^+2  sin  0  sin  0^008(0 )| 


(3. 15) 


for  kr  »  1,  &£•*,  and  0+0^  <  20q-t.  Higher  order  approximations  are  discussed 
by  Felsen  (1957d). 

i  -ikz 

When  a  plane  wave  is  incident  along  the  axis  of  symmetry,  i>  =  e  ,  the  total 
field  is  given  in  Section  n  (see  also  Siegel  and  Alperin,  1952  and  Felsen,  1955)  as 


2J>  -*»/2  l,(tr)^(c°ae) 

*  jW* 

"o 


(3. 16) 
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and  that  (3. 19)  and  (3. 20)  differ  only  In  sign  from  the  corresponding  soft  cone  re¬ 
sults  (3. 8  and  3.  9)  plotted  in  Fig.  5. 

Corresponding  to  (3. 17)  and  (3. 18),  the  following  expressions  for  the  back 
scattering  cross  section  obtain: 

X2(*-80)4  . 

«2«>-  "64,  (3  2D 

when  and 

n 

o2<0)= - - - -  (3.22) 

16*<0o-ir/2)4 

when  d0*  */2. 

A  more  general  expression  for  the  back  scattering  cross  section  is  given  by 
Schensted  (1953)  and  Siegel  et  ar(1953a,  1955b)  as 


where  the  summation  is  over  the  non-negative  zeros  of  (cos  6  ),  but  only  a  finite 

do  o 

o 

number  of  terms  must  be  included  since  the  infinite  series  diverges.  (See  remarks 
following  equation  3.  32).  Despite  this  drawback,  the  expression  was  used 
(Schensted,  1953;  Siegel  et  al  1953a,  1955b)  to  yield  the  next  order  in  the  large  and 
small  cone  angle  approximations  (3.  21)  and  (3. 22).  These  are 
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»2<°>-  -lXHco.90)] 


(3.24) 


for  6  ~  ic,  and 
o  ’ 


X2(1-4cos26  ) 
cr2(0)-  4 

16 t  cob  e0 


(3.25) 


for  9  ir/2. 

o  ' 


The  cross  sections  given  in  (3.  24)  and  (3.  25),  as  well  as  the  first  order  approxi¬ 
mations  (3. 21)  and  (3. 22),  are  plotted  in  Fig.  7  as  functions  of  60. 


3.  3  Vector  Boundary  Condition. 


The  field  of  a  transverse  electric  dipole  of  moment  4»  ejj^  and  arbitrary 


location  is  given  by 


fcfe-I 


I r-rl I 


(3.26) 


A  Q 

where  p^=p^  ®1+^  *1  •  When  this  source  is  in  the  presence  of  a  perfectly  con¬ 
ducting  semi-infinite  cone,  Felsen  (1957b)  gives  the  total  field  as 


I  J\  8  3i  % 

“  A/M  \.rl  80larl  +  9^9^ 


LA  *1  2\l 

n>,  ‘  VU10J  afjb 


(3.27) 


59 


THE  UNIVERSITY  OF  MICHIGAN 

3648 -2 -T 


6i 


THE  UNIVERSITY  OF  MICHIGAN 

3648 -2 -T 


The  notation  is  that  defined  under  equation  (3. 1).  Arbitrarily-oriented  elec¬ 
tric  dipoles  may  always  be  decomposed  into  transverse  and  radial  components,  so 
that  equation  (3. 27)  together  with  (3. 1)  is  sufficient  to  describe  the  field  exterior  to 
a  perfectly  conducting  cone  in  the  presence  of  an  electric  dipole,  regardless  of  its 
location  or  orientation.  An  alternative  representation,  in  terms  of  the  tensor 
Green's  function,  is  contained  in  Felsen  (1957a),  which  also  provides  a  comparable 
expression  for  an  arbitrary  magnetic  dipole  source. 

When  A  -ikr(sin  0  sin  0 

E  =  ^  e 

the  scattered  far  field  is  given  by  Felsen  (1957b)  as 


1  cos(0  -  0^ )  +  cos  0  cos  0^) 


Q 

ikr  /x-0  \  fsec^  §  +  see**  sin20  +  sin20. 

E0~  lk78in(l*~,4l)\  2/  cos 0+ cos  0  +  “  Td  2  77  ,2 1 

L  1  2 cos  -cos  -^(cos0  +  cosO^ J 


,  ikr  /jr-0  \ 

3=e_  ( - o) 

I  ikr  V  2  / 


2  sin0  sind 


(cos0  +COS0| ) 


sin20+sin20, 


g-  +  cos(0-jM 


r  2  0  2  01 

sec  —  +  sec  — 


cos  0  +  cos  0, 


16sin2  f sin2 


2cos2  |-cos2-tj(cos0+cos0^  (cos0+cos0i)2  J 


1 


(3.29) 


under  the  restrictions 

kr»  1,  0o">t,  0+01<20q  -j  . 
Note  that  ^  =  when  ^  =  -ir/2  . 
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to  the  particular  case  of  back  scattering  (9  =0^,  0  =^>  the  above  expressions 


for  the  field  components  reduce  to 

i  5  -ikr 
E  =^e 

E®=0 


ikr  /it -6  \  2. 

,  s  _  e _ / _ o  \  3+cos  8 

ikr  V  2  /  .  3 

r  >  '  4cos  0 


(3.  30) 


When  a  plane  wave  is  incident  along  the  axis  of  symmetry,  that  is. 


_i  a  -ikz 
E  =  i  e 
~  x 


(3.31) 


the  scattered  far  field  is  given  by  Hansen  and  Schiff  (1948)  (see  also  Mentzer,  1955) 


-lVpT 


~iv$*  .  -‘V  a  i 

.  e  P‘  (cob8)  e  j-l?v  (cos  6) 

s  ie^  *  V  &  66  a 

j0=._CO8^^  — - - 2_  ~-x - 

^  Binej  Q»^(/iJ2d^  va  ^  [p^(/^2 


.  ikr  r-1 

e5=  sin0  ■  L, 

V  kr  v. 


~ivf  a  1 
6  i?^{oose) 


-V  i 

e  P„  (cos© ) 


-Zs 


(3.  32) 


where  vQ  is  such  that  P*  (a»o)=0,  va  >4/2,  is  such  that  —■  1^(mo)=0,  i^>  -1/2, 
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and  kr  »  1 . 


No  upper  limit  is  given  for  the  sums.  This  omission  stems  from  the  fact 
that  the  vQ  and  are  infinite  sets  and  if  all  the  terms  are  included,  the  sums 
diverge.  Siegel  et  al  (1953a)  estimate  that  when  6  = 0  (back  scattering)  approximately 


"JLl 

_  100  J 


terms  are  of  the  form  shown  but  the  remaining  terms  (infinite  in  number) 


are  negligible.  This  complication  is  avoided  in  the  integral  representation  of  the 
scattered  field  given  by  Felsen  (1957b,  d).  Goryanov  (1961)  explicitly  gives  the  in¬ 
tegral  representations  corresponding  to  (3. 32)  as 


.  ikr 
8  ie 
E  =  ~  - —  < 
8  kr 


)-00  -00 

j  i/sinh(i/*)F1(i/)di/+-^^  j  v  sinh(va)F2(v)dv 


2  corf*  ^ 


(3.33) 


ilcjr  /*co 

8=ii__Bln0  J  „  sinh^^F^ ^  j  V  sinhMF^dv 

A  A  _  1 


2  °0 
cot2  2~ 

2  cos^  | 


where 


IW-V,2('CO80o)Piv^J<COBe) 

F>)=  -T - 

Piv  i/  i+  V4) 
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4Pl-1/,<’oo,8‘>)pi'-vl(co'9) 

*>w- 1 - 


ae  p  \v-^cobQo )cosh<,,*M*/2+ V4> 


and 


kr  »  1. 

Equation  (3.  33)  is  valid  only  for  0  <  20q-t,  that  is  for  field  points  exterior 
to  the  cone  formed  by  the  specularly  reflected  rays.  For  this  range  Goryanov 
(1961)  reports  extensive  calculations  of  the  bracketed  terms  in  (3. 33).  Although  his 
numerical  results  are  presented  only  qualitatively,  he  concludes  that  the  difference 
between  the  true  far  field  and  that  predicted  by  physical  optics  (see  below)  increases 
uniformly  with  angle  away  from  the  back  scattered  direction  and  that  the  difference 
never  exceeds  10  percent  of  the  value  computed  with  (3. 33)  for  0  <  20Q-ir— ar/  90. 

The  far  field  components  for  the  plane  wave  (3. 31)  incident  along  the  axis 
of  symmetry,  as  computed  by  physical  optics  (see  section  IV)  are  given  by 
Goryanov  in  the  form 


where 


ikr  ikr 

E  =  -  — — cos  f)  L,  E?=  — — sint)  L 

6  kr  r  J9  kr  r 

o 

sin  0ocos  dQ 

4cos  0/2|cos(0o-0/2)cos(0o+  0/2)]  ^ 


(3. 34) 


The  back  scattered  far  field  for  this  same  incident  field  (3. 31)  has  been 
approximated  by  Felsen  (1953, 1955)  as  follows: 
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ES( 0,0,2)=  -i  (l  -2 

X  V  *41-1)1 


(3.  35) 


where 


?  =  |  fkz  cosflj 
en'Jt/2>  kz»l,  C<4; 


(3.  36) 


0  ~  t/2,  kz  »  1,  ?  >  4; 


8  *  ie0'*  fw*°\ 


(3.  37) 


»  kr  »  1. 


Note  that  the  magnitudes  of  E  given  in  (3. 35)  and  (3.  36)  are  plotted  in  Fig.  5. 

Although  all  the  field  expressions  given  above  can  be  used  to  calculate 
radar  cross  sections,  the  calculations  which  have  been  carried  out  explicitly  in¬ 
variably  deal  with  plane  wave  Incidence  along  the  axis  of  symmetry.  The  physical 
optics  result  is  then  independent  of  f>  and  is  given  by  Siegel  et  al  (1955a)  as 


\2tan40 


a(0)= 


2(l+cos20  r 
_  o 


16t  (1+cos0Xcos0+cos20^3 


(3.  38) 
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When  6  =  0  this  reduces  to  the  widely-quoted  back  scattering  physical  optics  cross 
section,  given  first*  by  Spencer  (1951)  as 


0(0)= -^tanVflo)  • 


(3.  39) 


The  small  cone  angle  approximation  of  Felsen  (1955)  (equation  3.  37)  yields 
substantially  the  same  cross  section: 


«<°>-  iir(,,V 


(3.40) 


»  kr  »  1  . 


The  large  cone  angle  approximation  (3. 36)  also  gives  essentially  the  same 
result  as  physical  optics,  namely 


•«»•  "kt  «w«4 


(3.41) 


for  0  *  */2,  kr  »  1,  ^kr  cos  0Q »  1  . 

The  back  scattering  cross  section  obtained  from  the  field  expressions 
(3. 32)  is  given  by  Hansen  and  Schiff  (1948)  (see  also  Siegel  et  al,  1953a,  1955b; 
andMentzer,  1955)  as 


'*  *  1  PM '*• 


-to# 


ygyi)e 


\  a 


V  p 

o 


ansen  and  8chiff(1948)  quote  an  earlier  version  of  Spencer's  work,  now  unavailable. 
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1  9  1 

where  va  and  are  the  zeros  of  PvJ,nQ)  and  Pv^  (*i0)  respectively  and  the 


same  convergence  difficulties  mentioned  under  (3. 32)  are  still  present.  Despite 
these  difficulties  the  expression  was  successfully  employed  by  Schensted  (1953) 

(see  also  Siegel  et  al,  1953a,  1955b;  and  Felsen,  1955)  to  yield  second  order  cor¬ 
rections  to  equations  (3. 40)  and  (3. 41).  These  are 

a  (0)= - - - (1-2cos2®0)  (3. 43) 

16»cos*0o 


for  0  >vir/2,  and 
o 

2 

a (0)=  ^  (l+cos0o)2  [l+3(l+cos®0)j  (3. 44) 

for  6  it . 
o 

The  cross  sections  given  by  equations  (3. 39),  (3. 40),  (3. 41),  (3. 43)  and 
(3. 44)  are  plotted  in  Fig.  9  as  functions  of  0Q  . 
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1  iGURE  9:  NOSE-ON  BACK  SCATTERING  CROSS  SECTION  OF  A  PERFECTLY 
CONDUCTING  CONE 
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IV 

APPROXIMATE  RESULTS  FOR  A  FINITE  CONE 

Even  for  the  semi- Infinite  cone  the  complexity  of  the  exact  solution  is  such 
that  detailed  computations  are  extremely  difficult,  and  some  approximation  is  al¬ 
most  essential  if  the  physical  content  of  the  solution  is  to  be  apparent.  For  the 
finite  cone  the  situation  is  still  worse.  The  only  exact  solution  available  is  for  a 
spherically  capped  cone  (Northover,  1962;  Rogers  et  al,  1962)  and  the  results  are 
not  yet  amenable  to  computation.  It  is  therefore  appropriate  to  investigate  the  ex¬ 
tent  to  which  approximate  techniques  are  applicable  to  this  type  of  body. 

In  this  section  we  shall  be  concerned  exclusively  with  the  back  scattering 
cross  section  for  plane  wave  incidence  on  a  circular  cone  with  any  of  several  dif¬ 
ferent  types  of  terminations.  The  low  frequency  behavior  is  discussed  in  Section 
4. 1  and  in  the  next  section  the  physical  optics  method  is  used  to  derive  the  scattered 
Reid  at  high  frequencies.  The  intermediate  range  of  frequencies  is  again  the  diffi¬ 
cult  one,  but  in  Section  4. 3  some  of  the  methods  for  extending  the  high  frequency 
results  into  this  range  are  presented. 

4.1  Low  Frequency  Scattering 

Consider  a  plane  electromagnetic  wave  incident  upon  a  body  all  of  whose 
dimensions  are  finite.  At  sufficiently  low  frequencies,  the  complex  amplitude  of  the 
scattered  electric  field  in  the  far  zone  oan  be  expanded  in  a  series  of  increasing 
positive  powers  of  k,  where  k  is  the  wave  number.  If  this  amplitude  is  denoted 
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by  S,  and  is  defined  as  the  coefficient  of  e  /kr  where  (r,  6, ft)  are  spherical  polar 
coordinates  referred  to  an  origin  In  a  neighborhood  of  the  body,  then 

CD 

S  =  k3  51  knf  (4.1) 

^0  “n 

where  the  are  vector  functions  of  the  angular  variables  0  and  ft,  the  body  para¬ 
meters  and  the  characteristics  of  the  incident  field,  but  are  independent  of  k.  Since 
S  is  a  dimensionless  function,  it  is  clear  that  each  power  of  k  must  have  associated 
with  it  a  like  power  of  some  body  dimension./  ,  and  for  sufficiently  small  values  of 
k/  the  series  in  (4. 1)  is  convergent.  The  leading  term  corresponds  to  Rayleigh 
scattering  and  when  substituted  into  the  formula  for  the  crosB  section  (see  introduc¬ 
tion  to  Section  3),  we  have 

a  =  4jk4|fJ2  (4.2) 

which  depends  on  the  wavelength  through  the  inverse  fourth  power.  This  is  the  well 
known  Rayleigh  law  and  holds  for  all  angles  of  scattering  including  forward.  For 
practical  purposes,  however,  the  case  of  back  scattering  is  of  most  interest  and 
attention  will  be  confined  to  this  for  the  remainder  of  the  section. 

The  simplest  body  is  the  sphere  and  for  this  it  is  a  straightforward  matter 
to  determine  as  many  of  the  as  desired  by  appropriate  expansion  of  the  Mle 
series  (Goodrich  et  al,  1961).  In  particular,  the  leading  coefficient  £Q  is  made  up 
of  contributions  from  an  electric  and  a  magnetic  dipole  with  relative  magnitudes  1 
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2 

and  1/2  respectively,  giving  a  net  contribution  3a  / 2  where  a  is  the  radius  of  the 
sphere.  Hence 

a  -  9»k4a6  (4.3) 

which  can  be  written  as 

•  *  & k4y2  •  <4-4> 

where  V  is  the  volume,  and  we  remark  that  to  this  order  in  k  any  material  of  non¬ 
zero  conductivity  appears  perfectly  conducting. 

Since  exact  results  (such  as  the  Mie  series)  are  known  for  very  few  bodies  it 
is  fortunate  that  the  f^  can  be  obtained  in  an  alternative  manner  by  solving  certain 
scalar  problems  in  potential  theory.  The  technique  is  described  in  Senior  and  Dar¬ 
ling  (1963)  and  has  been  applied  to  the  spherically  capped  cone  by  Darling  (1960), 
leading  to  an  expression  for  1^  in  terms  of  infinite  matrices.  The  method  is  for¬ 
mally  applicable  to  any  body  which  is  the  intersection  of  a  finite  number  of  regions 
within  each  of  which  the  electrostatic  Green's  function  is  known,  and  can  therefore 
be  used  in  the  treatment  of  cones  with  other  types  of  termination. 

For  general  bodies  of  revolution  with  on-axis  excitation  an  approximate 
method  for  the  determination  of  |  f  |,  and  hence  of  o,  has  been  proposed  by  Siegel 
(1959).  The  theory  has  its  origin  in  the  fact  that  as  regards  fQ,  the  body  appears 
only  as  a  singularity  to  the  wave,  and  the  structure  of  the  body  is  in  the  nature  of  a 
second  order  correction.  Since  the  return  is  also  insensitive  to  changes  in  polari¬ 
zation,  fQ  can  be  a  function  only  of  the  geometric  properties  of  the  shape,  and  the 
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obvious  first  choice  for  a  function  of  dimension  J.  is  to  take  It  proportional  to  the 
volume.  Accordingly,  we  write 

o  =  C1k4(VF)2  ,  (4. 5) 

where  is  a  pure  number  and  F  is  a  dimensionless  correction  factor  that  takes 
into  account  the  approximate  dependence  of  the  cross  section  on  the  shape.  From  a 
study  of  the  exact  expression  for  |fQ  j  for  a  spheroid  of  semi-major  and  semi-minor 
axes  a  and  b  respectively,  Seigel  deduced  an  expression  for  F  of  the  form 

F  -  1  +  —  e"y  ,  (4.6) 

wy 

where  y  =  b/a,  and  with  =4/*  the  resulting  formula  for  o  is  in  agreement  with 
the  exact  spheroid  answer  (Rayleigh,  1897)  to  within  a  few  percent  for  all  elliptici- 
tles.  For  the  particular  case  of  the  sphere,  y  -  1  and  equations  (4. 5)  and  (4. 6)  then 
imply 

o  =  l.  5889  kSr2  , 

compared  with  the  correct  result  (equation  4. 4) 

o  =  1.6114  kSr2. 

If  it  is  now  assumed  that  (4. 5)  and  (4. 6)  are  a  valid  approximation  for  all 
bodies  of  revolution  under  symmetrical  excitation  with  y  taken  as  a  measure  of  the 
elongation  (the  ratio  of  a  characteristic  dimension  along  the  axis  to  one  in  the  per¬ 
pendicular  direction),  the  formula  can  be  used  to  derive  the  Rayleigh  scattering 
cross  section  of  a  finite  cone.  Inasmuch  as  the  nature  of  the  termination  is  unlikely 
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to  be  important  to  this  order  of  approximation*  it  ia  sufficient  to  oonfine  attention  to 
a  right  circular  cone  of  altitude  h  and  base  radius  a.  'Hie  volume  is  then  ra  h/3, 
and  from  the  definition  of  y  it  would  seem  natural  to  choose  y  *  h/2a.  A  better 
approximation,  however,  is  to  require  that  the  cross  section  approach  that  of  a  disc 
as  h— »  0,  from  which  we  have  y  *  h/4a  (Siegel,  1959).  The  cross  section  of  the 
cone  in  the  low  frequency  limit  is  therefore 

2 

4x  4  4.2  j  4ae-h/4a~l 

®  Tk  * b  I1  f  ,h  J  <4  ,) 

The  above  formula  is  in  good  agreement  with  experiment  (Brysk  et  al,  1959), 
though  care  is  necessary  to  ensure  that  the  frequency  is  low  enough  to  justify  its  use. 
Since  the  cross  section  is  determined  by  the  volume  and  the  maximum  dimensions, 
the  upper  extremity  of  the  Rayleigh  region  for  moBt  right  circular  cones  is  dictated 
by  the  slant  length  and  not  the  base  diameter  even  at  nose-on  incidence,  which  im¬ 
plies  that  for  a  cone  of  fixed  diameter  the  Rayleigh  region  recedes  with  decreasing 
cone  angle.  In  order  to  use  equation  (4. 7)  we  therefore  need 

k  Vh2+  r2  «  1 

and  is  the  reason  the  experimental  data  of  Keys  and  Primich  (1959a)  failed  to 
display  the  anticipated  wavelength  dependence.  We  also  observe  that  equations  (4. 5) 
and  (4. 6)  are  the  same  whether  we  view  the  object  from  the  front  or  the  rear,  and 

^Except  for  the  extent  to  which  the  volume  is  then  changed;  but  see  Siegel  (1963). 
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consequently  for  a  cone  in  the  Rayleigh  region  the  nose-on  and  base-on  cross  sec¬ 
tions  should  be  identical.  This  was  first  pointed  out  by  Siegel  et  al  (1959)  and  has 
since  been  verified  experimentally. 

So  far  we  have  considered  only  electromagnetic  scattering  but  a  few  words 

are  necessary  about  the  corresponding  scalar  problems  associated  with  hard  and 

soft  bodies.  The  hard  body  (Neumann  boundary  condition  at  the  surface)  is  the  one 

of  interest  in  acoustics  and  for  this  the  low  frequency  expansion  is  similar  to  that 

shown  in  equation  (4. 1)  except  that  the  coefficients  f  are  now  scalar  functions. 

n 

Only  a  single  dipole  contributes  to  fQ,  however,  and  the  exact  Rayleigh  cross  sec¬ 
tion  for  a  hard  sphere  is 


(4.8) 


which  is  smaller  that  the  electromagnetic  result  by  a  factor  81/25.  Siegel's  approxi¬ 
mate  method  leads  to  a  value  =  l/w  (a  single  dipole  effect  rather  than  two  equal 
in-phase  contributions)  so  that  for  a  hard  cone  the  cross  section  given  in  (4. 7)  must 
be  reduced  by  a  factor  4.  For  a  soft  body,  on  the  other  hand,  the  expansion  for  the 
far  field  amplitude  S  is  quite  different.  The  leading  term  is  proportional  to  k  with 
a  coefficient  which  is  simply  the  electrostatic  capacity  of  the  body  (Senior  and  Dar¬ 
ling,  1963),  and  the  entire  scattering  behavior  now  differs  accordingly. 

4.2  High  Frequwnny  Spattering 

At  high  frequencies  the  scattering  behavior  is  more  intimately  related  to  the 
(geometrical  properties  of  the  body,  and  if  the  wavelength  is  small  in  comparison 
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with  all  the  major  dimensions,  an  accurate  approximation  to  the  scattering  cross 
section  can  be  obtained  by  the  method  of  physical  optics.  In  many  cases  the  method 
gives  good  results  even  when  there  are  singularities  cm  the  surface,  such  as  edges 
and  points  (as,  for  example,  the  vertex  of  the  cone),  and  we  shall  begin  by  using  it 

i  "Iky,  a 

to  calculate  the  back  scattering  cross  section  for  a  plane  wave  (E  =  e  i  )  at 
nose-on  incidence  on  a  cone  with  various  types  of  termination. 

The  method  is  based  on  an  approximation  to  the  field  distribution  induced  on 
the  surface  of  the  object.  On  the  illuminated  side  of  the  surface  each  element  is 
assumed  to  possess  the  field  it  would  have  possessed  had  it  been  part  of  an  infinite 
tangent  plane,  whereas  on  the  shadowed  portion  the  field  is  taken  to  be  zero.  This 
leads  to  a  current  distribution  which  is  zero  in  the  shadow  and  twice  the  tangential 
component  of  the  incident  magnetic  field  in  the  lit  region,  and  the  calculation  of  the 
scattered  field  is  now  reduced  to  quadratures. 

The  general  expression  for  the  scattered  field  under  this  approximation  is 
given  in  Crispin  et  al  (1959).  This  simplifies  considerably  for  the  case  of  back  scat¬ 
tering  from  a  body  of  revolution  at  nose-on  incidence.  If  such  a  surface  is  defined 
by  the  equation 

z  =  g(p)  ,  p  =l|x2  +  y2 

then  the  scattered  electric  vector  in  the  far  field  can  be  written  as 
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E8  *  ESf  , 
-  x 


ES=  * 


ikr 


kr 


^ik2  ^  e~2ikgpdp) 


(4.9) 


where  the  Integration  is  carried  out  over  the  illuminated  portion  of  the  body  (for  a 
smooth  convex  body  a  =  0,  b  =  Pg^dow^'  ^he  bracketed  quantity  in  (4. 9)  is  simply 
the  far  field  amplitude  and  if  this  is  denoted  by  S  then 


c  4x2i  \  -2ikg  . 

S  =  — —  \  e  pdp 

*  Ja 


in  terms  of  which  the  cross  section  is  given  by 


For  a  flat-backed  cone  of  height  h  and  semi-vertex  angle  r  - 


c  41t  2a  \  2iku  . 

S  =  — —  tan  6  \  e  udu 

X2  0 


(4.10) 


(4.11) 


=  -■>  tan2eQ  jl  +  e2Uch  (2ikh  -  1)J-  , 


(4.12) 


which  implies  a  scattering  amplitude 


tan20 
4  o 


(4.13) 
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associated  with  the  tip,  and  a  scattering  amplitude 


1.2  .  /011 .  2ikh 
-  —  tan  9  (21kh  -  l)e 
4  0 


(4. 14) 


corresponding  to  the  edge  singularity  at  the  base.  For  large  (but  finite)  kh  the 
latter  is  approximately 


xh ,  2  2ikh 
—  tan  9  e 
X  o 


(4.15) 


which  dominates  the  overall  return  at  high  frequencies  due  to  its  wavelength  depen¬ 
dence.  In  practice,  however,  a  more  accurate  estimate  of  the  base  return  is  re¬ 
quired,  and  this  is  considered  later. 

In  the  limit  of  a  semi-infinite  cone  (h  =  oo)  no  return  can  come  from  the  "ter¬ 
mination"  and  the  second  of  the  above  contributions  muBt  be  ignored.  This  is  usually 
justified  on  a  mathematical  basis  by  attributing  a  non-zero  conductivity  to  the  sur¬ 
rounding  medium  (giving  k  a  positive  imaginary  part),  so  that  the  second  term  in 
(4.12)  is  exponentially  attenuated  through  the  factor  e2^*.  The  only  return  is  now 
provided  by  the  tip,  and  from  (4. 13)  we  have 


X  ,  4- 
a  *  —  tan  9 

IDT  O 


(4.16) 


which  is  the  well  known  "tip"  answer. 

With  terminations  other  than  a  flat  back  the  physical  optics  result  will  only  be 
affected  if  a  portion  of  the  cap  is  in  the  illuminated  region,  and  consequently  for  a 
spherically-capped  cone  with  the  origin  of  the  sphere  at  the  apex  of  the  cone  (or. 
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Indeed,  on  the  aids  of  the  cone  at  any  distance  X  from  the  apex  with  0  <  X  <  h)  the 
predicted  return  is  identical  to  that  in  equation  (4.12).  But  if  the  sphere  is  chosen 
to  have  its  origin  at  a  distance  h  +  X  (X  >  0)  from  the  apex,  with  radius 
■{(htan  0q)‘  \jyi2  (so  as  to  achieve  a  Join  to  the  cone  with  a  discontinuity  in  at 
most  the  first  derivative  of  the  profile*),  a  part  of  the  sphere  is  visible,  and  the 
physical  optics  cross  section  is  changed  accordingly.  The  far  field  amplitude  is  now 

S  =  -  ~  jtan20Q  +  e2ikh  j2ik(htan20Q-  X)  -  sec^J  +  e2ik(x+h)J.  (4. 17) 

The  first  term  arises  from  the  tip  and  is  as  before.  The  second  group  of  terms 
represents  the  scattering  amplitude  associated  with  the  join,  and  can  be  written 
approximately  as 

2ikh 


S  =  7  sec2© 
4  o 


(4. 18) 


jl  -  2ikhtan  0Qtan<S  j 

where  S  is  the  angle  between  the  tangents  to  the  cone  and  the  sphere  at  their  Junc¬ 
tion.  Hence,  for  the  cone- sphere  (6  =  0)  the  scattering  amplitude  attributable  to 
the  join  is 

<  O  041*4% 

(4.19) 


0  i  2  21kh 
S  *  7  sec  0  e 
4  o 


which  has  the  same  wavelength  dependence  as  the  tip  contribution,  but  is  of  a  much 
larger  magnitude  for  all  cones  except  those  of  largest  angle.  It  is  believed  that  this 


The  profile  is  simply  the  curve  z  =  g (p)  which  defines  the  body. 
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is  an  accurate  estimate  of  the  return  when  the  first  derivatives*  of  the  profile  are 
matched.  Inasmuch  as  the  surface  in  the  neighborhood  of  the  singularity  is  entirely 
in  the  illuminated  region,  the  assumed  current  distribution  should  not  be  appreciably 
in  error,  and  the  most  accurate  experimental  data  available  in  the  Radiation  Labora¬ 
tory  is  consistent  with  a  scattering  amplitude  having  this  wavelength  dependence, 
and  a  magnitude  differing  from  that  in  (4. 19)  by  less  than  20  percent. 

On  the  other  hand,  if  the  kink  angle  5  is  not  zero,  the  second  term  in  (4. 18) 
is  present,  and  will  dominate  the  return  at  high  frequencies.  The  corresponding  far 
field  amplitude  is 

8  *  sec^0  tanS  e^*^  ,  (4.20) 

X  o 

where  a  is  the  radius  of  the  singularity,  and  this  is  almost  identical  to  the  result 
obtained  by  Dawson  et  al  (1960)  using  circular  wedge  theory  (see  later  discussion). 

It  is  therefore  obvious  that  in  any  experimental  measurement  of  cone-sphere  cross 
sections  extreme  care  must  be  taken  to  ensure  accurate  modeling  of  the  surface  pro¬ 
file  near  to  the  cone-sphere  Join. 

The  third  term  in  equation  (4. 17)  originates  at  the  shadow  boundary  of  the 
sphere  and  is  known  to  be  spurious.  It  arises  because  of  the  discontinuity  in  the  cur¬ 
rent  distribution  introduced  by  the  physical  optics  approximation,  and  is  independent 

*The  wavelength  dependence  of  the  physical  optics  contribution  from  a  Junction  be¬ 
tween  two  surfaces  is  Xn“^,  where  n  is  the  order  of  the  first  discontinuous  deriva¬ 
tive. 

J 
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of  the  presence  of  the  cone.  From  an  analysis  of  the  field  scattered  by  a  sphere 
(see,  for  example,  Goodrich  et  al,  1961),  It  Is  seen  that  the  only  contribution  from 
the  shadow  boundary  Is  produced  by  a  set  of  creeping  waves  which  have  traversed 
the  back  of  the  sphere.  The  associated  scattering  amplitude  decreases  exponentially 
with  increasing  frequency,  and  is  small  compared  with  the  specular  return  for  all 
spheres  more  than  a  few  wavelengths  in  radius.  With  the  cone-sphere,  however, 
the  specular  return  is  not  present,  and  since  the  contributions  from  the  tip  and  join 
are  quite  small  in  magnitude  (and  certainly  so  when  compared  with  a  specular  signal) 
the  creeping  waves  continue  to  provide  a  significant  effect  until  the  wavelength  has 
decreased  to  a  very  small  fraction  of  the  sphere  radius.  Nevertheless,  at  suffi¬ 
ciently  high  frequencies  the  creeping  waves  can  be  ignored,  and  there  is  some  slight 
evidence  (Turner  and  Dawson+,  1960;  unpublished  Radiation  Laboratory  data)  to  sug¬ 
gest  that  the  best  estimate  of  the  cone-sphere  return  is  then  obtained  by  ignoring  the 
spurious  term  in  (4.17).  The  resulting  high  frequency  cross  section  is 

4  . 

o  =  —t—  (tan  Q  +  sec  6  )  (4.21) 

lox  o  o 

and  since  the  term  produced  by  the  phase  factor  e^^  has  been  omitted,  (4.21) 

should  be  interpreted  as  an  average  over  a  small  frequency  band.  For  small  cone 

2 

angles  the  cross  section  is  of  order  0. 02  X  ,  which  is  approximately  one-half  that 
which  would  have  been  found  if  the  spurious  term  had  been  retained 

+who  also  use  the  physical  optics  method  to  calculate  the  return  due  to  a  termination 
produced  by  the  rotation  of  an  off-set  spherical  arc  about  the  cone  axis. 
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Before  leaving  this  discussion  of  the  nose-on  behavior,  a  final  point  of  inter¬ 
est  is  the  transition  of  the  cross  section  for  the  flat-backed  cone  to  the  diso  answer. 
As  the  half-cone  angle  r  -  Qq  approaches  w/2,  with  the  base  radius  a  held  constant, 

the  cone  degenerates  into  a  disc  whose  far  field  amplitude  is  inversely  proportional 
2 

to  X  .  The  manner  In  which  this  transition  occurs  can  be  seen  from  equation  (4. 12) 
if  tan  is  replaced  by  -a/h.  We  then  have 

S  =  -  ~  (a/h)2  jl  +  e2ikh(2ikh  -  l)j  (4. 22) 

and  by  expanding  the  exponential  function  for  small  argument  it  is  found  that  as 
h->  0  the  bracketed  terms  decrease  in  just  such  a  way  as  to  cancel  the  infinity  pro¬ 
duced  by  the  first  factor.  In  the  limit  h  =  0  the  far  field  amplitude  is  simply 

S  =  2i(y^  (4.23) 

which  is  the  known  result  for  a  disc  of  radius  a.  We  note  in  passing  that  this  is  also 
the  appropriate  expression  for  rear-on  incidence  on  a  flat-backed  cone. 

For  angles  of  incidence  away  from  nose-on  the  calculation  of  the  back  scat¬ 
tering  cross  section  by  the  physical  optics  method  is  again  straightforward,  though 
more  tedious  owing  to  the  asymmetrical  illumination  and  (possibly)  the  varying  posi¬ 
tion  of  the  shadow  boundary  as  a  function  of  aspect.  Moreover,  because  of  the  inac¬ 
curacy  in  the  prediction  of  the  shadow  boundary  effect,  the  usefulness  of  the  method 
is  mainly  limited  to  the  determination  of  the  specular  return  from  the  cone  sides  and 
its  side  lobes.  For  such  a  calculation  it  is  generally  sufficient  to  assume  that  the 
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cone  has  a  flat  back,  though  the  type  of  termination  will  affect  the  manner  in  which 
the  return  falls  off  at  aspects  beyond  the  specular  direction. 

If  the  plane  wave  is  Incident  at  an  angle  a  to  the  negative  z  axis,  the  optics 
integral  can  be  obtained  from  the  general  expression  in  section  4. 3  of  Crispin  et  al 
(1959).  When  this  is  particularized  to  the  problem  of  a  right  circular  cone  and  the 
notation  changed  to  that  presently  in  use,  the  expression  for  the  far  field  amplitude 
becomes 


~^2~  ^*0  ^>inoreos0  -  cost* tan 0^ 


2ikz'(cosa  +  sinatan0  cob0)  ,, 

•  e  o  r  z'dz’d^ 


(4.24) 

where  z'  is  measured  along  the  axis  of  the  body  and  the  integration  is  over  the  illu¬ 
minated  portion.  The  result  is  independent  of  whether  the  electric  vector  lies  in  the 
plane  containing  the  axis  and  the  direction  of  propagation  (horizontal  polarization)  or 
in  the  plane  perpendicular  thereto  (vertical  polarization). 

When  the  direction  of  incidence  lies  within  the  backward  cone,  i.  e. 

0  <a  <  r  -  0  ,  the  whole  of  the  surface  is  illuminated  and  the  z'  integration  can  be 
o 


carried  out  immediately  to  give 


-r—  tan  0 
8*  o 


sina  cos  0  -  c  os  a  tan  0 
_  o 

2 

(cos  or  +  sina  tan  0  cos  f) 

1  ~T  0 


[-  0- 


2ikh(cosa 


+  sina  tan  0qcos  0)^  e' 


2ikh(cosa+sinatan0  cos  f>) 


(4.25) 
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In  the  special  case  of  nose-on  incidence  (or  *  0)  the  f)  integration  is  trivial,  leading 
to  the  result  shown  in  equation  (4. 12),  but  otherwise  the  validity  of  (4.25)  is  limited 
to  angles  a  away  from  the  specular  direction  t/2  (which  lies  within  the  back¬ 
ward  cone  only  for  half  angles  x  -  9q  >  t  /4).  The  integration  of  the  first  term  in 
(4. 25)  gives+ 


*  2a 
t  tan  6 

*  (1  -  sin2osec20  )3^2 


(4.26) 


which  is  the  required  generalization  of  the  tip  return  (4. 13).  If  in  addition 
khslnatan(x-  6q)  « 1,  the  second  group  of  terms  in  (4.25)  can  be  treated  by  ex¬ 
panding  the  exponential  for  small  argument,  leading  to  a  simple  multiple  of  (4. 26) 
but  with  a  phase  corresponding  to  a  contribution  from  the  base  of  the  cone.  Hence, 

for  incidence  in  the  backward  cone  with  khsinatan(x-  6  )  «  1, 

o 


S  “  -  -7  tan2  9  (l+3sin2asec20  ) 
4  o  0 


jjl  +e 


2ikhcosa 


(21khcosa 


"»] 


(4.27) 


which  should  be  compared  with  the  expression  for  nose-on  incidence  shown  in  equa¬ 
tion  (4. 12).  On  the  other  hand,  if  khslnatan(z  -  9  )  » 1  the  integration  of  the 

o 

second  group  of  terms  in  (4. 25)  can  be  carried  out  by  the  stationary  phase  method, 


and  the  resulting  far  field  amplitude  is 


I 

4 


{ 


o  .  ,  „  1  21khcos(ar- 6  )sec 8 

- : -  tan  (a  -  9  )  e  0  0 

xsinar  o 


A  result  first  obtained  by  Spenoer  (1951). 


(4.28) 
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The  phase  Is  now  appropriate  to  a  contribution  from  the  base  end  of  the  nearest 
generator  of  the  cone,  and  because  of  the  restriction  on  kh  under  which  the  above 
expression  was  obtained,  the  tip  contribution  is  negligible  by  comparison. 

For  angles  of  incidence  outside  the  backward  cone  only  a  portion  of  the  sur¬ 
face  is  illuminated,  but  this  fact  is  of  no  concern  if  the  stationary  phase  method  is 
applied  to  the  0  Integral  in  (4.24).  We  then  have 


{ 


tanOr-  0  ) 
o 

2Xslnor 


sin(or  -  0  )  sec  0 
o  o 


2ikz'coe<«  -  0  )sec  6 
a  °  °{? 


dz' 


which  can  be  written  as 


ie**/4 


i 


khtan(v-  6  ) 
o 

v  s  In  a 


2ikhcos(or-0  )sec0 


tan  (a  -  0  )  e 
0 


where 


{i  -  P  (v 


2khcos (pt -6  )sec  0 
o  o 


)} 


(4.29) 


F(t)  =  - 


-tr 


ir2  . 
e  dr , 


yJo 

and  is  related  to  the  Fresnel  integral  (see,  for  example,  Jahnke  and  Emde,  1945). 

If  1 1 1»  1,  F(C)  *  CXt'1)  and  can  be  neglected  in  comparison  with  the  first  term 
in  (4. 29),  in  which  case  we  recover  the  expression  shown  in  (4. 28).  For  small  |c| , 
however. 
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F(t)  =  1  -||it2  +  0(t4) 
and  consequently  as  a  approaches  t  / 2,  the  far  field  amplitude  reduces  to 


iv  *  t  m  I  khsinfl 

?'  7  l/T-2 


2 

sec  0  . 
o 


(4.30) 


This  is  the  specular  contribution.  It  is  equivalent  to  the  broadside  return  from  a 

thick  cylinder  whose  length  is  the  slant  length  of  the  cone  and  whose  radius  is 

4/9  asec(T  -  9  ),  where  a  =  htan(r  -  6  )  is  the  base  radius  of  the  cone,  but  this 
o  o 

choice  of  dimension  is  not  unique,  nor  is  it  sufficient  for  computing  the  side  lobes. 
Indeed,  the  far  side  lobes  must  be  calculated  from  equation  (4.28)  and  to  investigate 
the  transition  from  (4.28)  to  (4.30)  it  is  necessary  to  determine  S  from  (4.29)  using 
the  available  tabulations  of  the  Fresnel  integral.  The  asymmetrical  nature  of  the 
side  lobes  is  then  apparent  even  for  the  flat-backed  cone,  and  can  be  observed  in 

the  experimental  data  of  Keys  and  Primich  (1959c,  d). 

As  previously  remarked,  one  of  the  limitations  of  the  physical  optics  approach 
is  the  inaccuracy  of  its  prediction  of  the  return  from  a  shadow  boundary.  If  a  flat- 
backed  cone  (or  any  finite  cone  whose  termination  leads  to  a  "positive"  discontinuity 
in  the  first  derivative  of  the  surface  profile  at  the  join)  is  viewed  at  or  near  nose-on 
the  shadow  boundary  coincides  with  the  base  of  the  cone,  and  the  physical  optics 
result  for  the  scattering  from  the  base  is  suspect.  Moreover,  at  high  frequencies  the 


This  is  also  evident  from  equation  (4.29)  if  the  Fresnel  integral  is  approximated 
[for  small  values  of  its  argument. 
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base  has  the  appearance  of  a  ring  singularity  at  these  aspects,  suggesting  that  the 
corresponding  contribution  will  be  dominant;  even  physical  optics  displays  this  fact 
(see  equation  4. 15),  but  before  any  reliance  is  placed  upon  the  magnitude  which  it 
predicts  an  alternative  method  of  derivation  is  desirable.  This  is  particularly  so 
since  (4.15)  can  be  written  as 

21kacot(*  -  6  ) 

V  Q  A 

s  =  ■—  tan(r  -  G  )  e  ;  (4.31) 

A.  O 

for  fixed  a,  the  far  field  amplitude  tends  to  zero  with  decreasing  cone  angle,  where¬ 
as  the  analogy  with  the  ring  would  suggest  a  limit  which  is  finite  but  non-zero. 

To  provide  a  more  accurate  estimate  of  the  base  return  at  high  frequencies, 
a  method  was  proposed  by  Siegel  in  1957  (see  Siegel,  1959)  which  is  similar  to  one 
employed  by  Artmann  (1950)  in  his  analysis  of  the  thick  half  plane.  Taking,  for 
example,  a  flat-backed  cone  at  nose-on  incidence,  the  idea  is  to  treat  every  element 
of  the  edge  formed  by  the  sides  of  the  cone  and  the  base  as  a  portion  of  an  infinite 
wedge  on  which  the  field  is  incident  at  an  appropriate  angle.  By  using  the  known  sol¬ 
ution  for  such  a  wedge,  and  regarding  every  element  as  scattering  independently  of 
the  rest,  the  net  return  can  be  obtained  by  simple  integration  providing  due  account 
is  taken  of  the  varying  polarization  for  the  different  elements.  The  solution  origi¬ 
nally  developed  was  in  error  because  of  an  incorrect  Interpretation  of  the  polariza¬ 
tion  effect.  The  error  was  pointed  out  by  Keller  (1959)  and  duly  corrected  (Siegel 
et  al,  1959).  Simultaneously,  Keller  (1959;  but  see  1960)  published  an  alternative 
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derivation  based  on  his  geometrical  theory  of  diffraction  (Keller,  1957),  and  pre¬ 
sented  the  solutions*  for  all  aspects,  hard  and  soft  soalar  cases  as  well  as  vector. 
For  the  electromagnetic  problem  with  nose-on  incidence  the  far  field  amplitude  asso¬ 
ciated  with  the  base  return  is 


2ikacot(x  -  9  ) 
2 x  o 

cosec  —  e 
n 


(4.32) 


where 


11  *2  ~  T  ’  (433) 

and  compared  with  the  physical  optics  solution  (4.31)  which  tends  to  zero  as  9  — >x 

o 

(cone  of  vanishing  angle),  (4. 32)  is  non-zero  at  this  limit. 

A  graphical  comparison  of  the  two  solutions  aB  a  function  of  cone  angle  is 

given  in  Siegel  et  al  (1959)  and  clearly  shows  the  increasing  discrepancy  as  9  — >  v . 

o 

For  small  cone  angles,  the  previous  estimates  of  the  base  return  must  therefore  be 
abandoned  in  favor  of  that  obtained  by  wedge  theory.  Accordingly,  (4. 15)  and  (4. 31) 
must  be  replaced  by  (4. 32),  but  for  a  termination  which  produces  only  a  slight  kink 
in  the  profile  (or  no  kink  at  all),  wedge  theory  and  physical  optics  are  in  agreement, 
and  consequently  the  results  for  the  cone-sphere  are  unaffected. 


Two  errors  in  Keller  (1959,  1960)  were  noted  in  1960  (see  Crispin  et  al,  1963)  and 
corrected  by  Keller  (1961).  The  second  term  in  equation  (29)  of  Keller's  papers 
should  be  (cos  x/n  -  cos  3x/n)“  ,  and  the  legends  attached  to  the  two  curves  in 
Fig.  4  should  be  interchanged. 
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As  >  t/2  (the  disc),  both  (4. 31)  and  (4. 32)  tend  to  infinity,  and  do  so  at 
the  same  rate.  Neither  of  them  now  provides  a  satisfactory  estimate  of  the  base  re¬ 
turn.  On  the  other  hand,  we  have  already  seen  that  the  three-term  optics  approxi¬ 
mation  (4. 12)  is  successful  in  reproducing  the  correct  transitional  behavior,  and  for 
wide  angle  cones  (4. 12)  should  be  used  in  place  of  either  (4. 31)  or  (4. 32).  As  of  the 
moment  one  must  arbitrarily  choose  the  cone  angle  at  which  to  switch  from  (4. 32)  to 
(4. 12),  but  it  has  been  suggested  (Crispin  et  al,  1963)  that  if  h  »a,  (4.32)  should 
be  used;  if  kh  « 1  «  ka,  (4. 12);  and  if  h  and  a  are  comparable,  the  mean  of  the 
two  results. 

In  an  attempt  to  predict  the  high  frequency  scattering  with  an  accuracy  grea¬ 
ter  than  that  which  can  be  achieved  by  physical  optics,  the  cylindrical  current 
method  has  been  pioposed  by  Dawson  and  Turner  (1960).  This  also  is  based  on  an 
approximation  to  the  current  distribution  excited  on  the  surface,  and  as  such  is  simi¬ 
lar  to  physical  optics,  but  whereas  the  latter  assumes  that  each  element  scatters  as 
though  part  of  an  infinite  tangent  plane,  the  new  approach  postulates  a  current  which 
is  the  same  as  for  an  infinite  cylinder  having  axis  and  radius  identical  to  those  of  the 
body  at  the  point  in  question.  It  is  therefore  limited  to  bodies  of  revolution,  but  for¬ 
mally  at  least  it  is  "valid"  for  a  much  wider  range  of  radius  to  wavelength  ratios 
providing  that  the  cylinder  solution  is  used  in  its  exact  form.  If  the  incident  electric 
vector  is  normal  to  the  plane  containing  the  direction  of  propagation  and  the  axis  of 
the  body  (vertical  polarization),  the  far  field  amplitude  obtained  by  this  method  is 


THE  UNIVERSITY  OF  MICHIGAN 

3648-2-T 


q  2i  \  e  2ikz'cosa  ,  , 

8  ■  x  \  si  •  **' 


where 


oo 


J'  (kR  slnar) 
d  n 


-oo 


(lP 

H'  '  (kR  slna) 
n 


(4.34) 


(4. 35) 


and  R  is  the  radius  of  the  body  at  a  distance  z*  along  its  axis.  For  horizontal 
polarization  the  solution  differs  only  in  having  replaced  by  S^,  the  correspondlngj 
sum  over  the  undifferentiated  functions. 

With  all  except  the  most  trivial  bodies,  an  exact  analytical  evaluation  of  the 
above  integral  is  impossible,  and  a  digital  computer  is  almost  essential  if  the  full 
advantages  of  the  method  are  to  be  realized.  Nevertheless,  if  kRsina  » 1,  -S^ 
and  Sg  can  both  be  approximated  by  the  function 

\  el<»/4-2Ut.i n«> 

and  when  expressed  in  terms  of  the  electric  vector  the  scattering  amplitudes  for  ver-j 
tlcal  and  horizontal  polarization  are  now  identical.  This  is  the  same  type  of  polari¬ 
zation  independence  displayed  by  the  physical  optics  method.  The  resulting  expres¬ 
sion  for  S  is 

e‘*/4  Jfif  (4.36) 

with  R  =  R(z'),  and  this  has  been  evaluated  by  Dawson  and  Turner  (1960),  Turner 
(1960)  for  bodies  of  various  shapes.  In  particular,  for  the  flat-backed  cone 
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R  *  -z'  tan  8  ,  0  <  z'  <  h,  and  the  far  field  amplitude  associated  with  the  return 
o 

from  the  cone  sides  is  then  of  the  same  form  as  that  shown  in  equation  (4.29),  but  is 

smaller+  by  a  factor  (1  -  tan  6  cot  a).  This  implies  a  reduction  in  the  specular  con- 

o 

2 

trlbution  by  a  factor  sec  8  .  Analogous  results  for  the  cone-sphere  have  also  been 

o 

given  by  Dawson  and  Turner  (1960),  and  these  reveal  a  similar  discrepancy  with  the 
physical  optics  predictions. 

In  spite  of  the  functional  differences  between  the  two  sets  of  formulae,  the 
differences  in  their  numerical  values  are  only  significant  for  wide  angle  cones,  and 
the  experimental  data  available  at  the  moment  is  not  sufficient  to  show  which  of  the 
two  methods  is  more  accurate.  For  a  body  of  revolution  the  cylindrical  current 
method  would  seem  to  have  the  greater  potential,  but  is  possible  that  the  intrinsic 
accuracy  is  lost  in  the  approximation  to  and  S^.  Certainly  this  approximation 
prevents  any  consideration  of  nose-on  incidence,  and  for  this  range  of  aspects  we 
are  forced  to  rely  on  physical  optics,  or  the  other  high  frequency  approximations 
described  in  this  section. 


For  a  body  of  revolution  whose  dimensions  are  greater  than  a  wavelength  or  so  and 
whose  surface  is  half  in  shadow,  Turner  (1960)  claims  that  the  two  methods  are,  in 
fact,  equivalent.  However,  the  proof  appears  to  assume  that  R  is  independent  of  z', 
and  since  the  body  is  then  a  cylinder  of  constant  radius,  this  is  the  one  case  for 
which  the  conclusion  is  obvious. 
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4.3  The  Resonance  Region 

The  resonance  region  can  be  regarded  as  the  frequency  range  for  which  the 
wavelength  is  comparable  to  a  major  dimension  of  the  body.  It  is  intermediate  to  the 
ranges  for  which  the  Rayleigh  and  optical  (high  frequency)  laws  apply,  and  is  usually 
characterized  by  the  presence  of  large  oscillations  in  the  back  scattering  cross  sec¬ 
tion  as  a  function  of  wavelength.  It  also  poses  one  of  the  severest  challenges  to  the 
theoretician  since  it  is  here  the  the  transition  between  two  contrasting  types  of  scat¬ 
tering  behavior  takes  place. 

In  attempting  to  investigate  this  region  it  is  natural  to  consider  the  possibility 
of  extending  the  techniques  appropriate  to  the  neighboring  frequency  ranges.  By 
including  more  terms  in  the  low  frequency  expansion  (4. 1)  it  may  bo  feasible  to  pene¬ 
trate  the  region  to  some  small  extent,  but  since  the  series  has  a  finite  radius  of  con¬ 
vergence  (Senior  and  Darling,  1963)  this  procedure  is  strictly  limited  both  in  theory 
and  in  practice.  On  the  other  hand,  no  formal  restriction  of  this  type  exists  at  the 
high  frequency  end,  and  at  least  a  partial  coverage  of  the  resonance  region  can  be 
achieved  by  applying  the  pseudo  asymptotic  techniques  available  at  high  frequencies. 

In  the  case  of  a  flat-backed  cone  the  necessary  extension  can  be  obtained 
using  the  geometrical  theory  of  diffraction,  and  for  nose-on  incidence  the  appropri¬ 
ate  formulae  have  been  derived  by  Keller  (1960).  Whereas  the  wedge  solution  (4.32) 
corresponds  to  those  rays  which  are  scattered  by  the  base  and  return  directly  to  the 
receiver,  the  seoond  order  term  is  produced  by  rays  which  traverse  the  base  of  the 
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cone,  and  are  therefore  diffracted  twice  before  being  observed.  The  far  field  ampli¬ 
tude  correct  to  this  order  is  now 


S  -  -~ 


»a  2t 

— -  cosec —  e 
nX  n 


21kaoot(r-  6 


>  J  ,  ,  1  .  w 

■<  1  +-  sln- 
^  n  n 


21ka-  lx/4 


(*ka) 
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t  3t 

COB-  -  cos - 

n 


(cos~ 


00ST  1 

3*2  f 
cos-l  J 


(4.37) 


where  n  Is  defined  In  equation  (4. 33),  and  the  curves  shown  in  Figs.  10  through  20 
are  based  on  this  formula.  The  agreement  between  theory  and  experiment  is 
remarkably  good  even  down  to  values  of  ka  of  order  unity.  Since  the  exponential 
within  the  braces  in  (4. 37)  has  a  negative  coefficient  for  all  cone  angles,  the  maxima 
in  the  cross  section  occur  at 

ka  =  (m  +  5/8 )x  , 
m  *  0, 1, 2. . . ,  and  the  minima  when 

ka  *  (m  +  9/8 )r 

(see,  for  example.  Fig.  4  of  Keller,  1960).  The  curve  Intersects  that  of  the  Ray¬ 
leigh  cross  section  at  around  ka  =  1,  and  between  them  these  two  formulae  are 
sufficient  to  indicate  the  variation  of  the  nose-on  return  throughout  the  resonance 
region. 

For  angles  of  lncldenoe  away  from  nose-on  an  accurate  determination  of  the 
back  scattering  cross  section  of  a  flat-backed  oone  is  not  so  Important  for  most 
purposes,  and  because  of  this  the  formulae  have  received  less  attention.  In  particu¬ 
lar  the  geometrical  theory  of  diffraction  has  only  been  taken  as  far  as  the  first  order 
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(high  frequency)  term,  but  It  may  be  appropriate  to  summarize  the  formulae  and 
give  a  few  remarks  about  their  applicability. 

If  the  Incidence  is  at  an  angle  a  £  0  but  within  the  backward  cone  the  far 
field  amplitude  corresponding  to  (4. 32)  is 

S  *  u(a)  +  u(-a)  (4.38) 

(Keller,  1960)  where+ 


“(or)  -  "COt^\/ 


a 


Xsina 


lsr/4  -  2ka cosec  6  cos(or -9  )> 
*-  o  o  J 


(4.39) 


with  the  upper  sign  for  vertical  polarization  (or  soft  body)  and  the  lower  for  hori¬ 
zontal  polarization  (hard  body).  We  observe  that  u(a)— »a>  as  or— >0:  this  is  a  con¬ 
sequence  of  the  caustic  which  exists  along  the  continuation  of  the  cone's  axis. 

The  above  expression  also  holds  in  the  aspect  range  t/2<o<^  providing 
3f  +  2a  2a  —  f 

the  term  cos -  is  replaced  by  cos - ,  but  is  infinite  at  a  ~  j .  This  is 

n  n 

again  attributable  to  a  caustic  and  to  the  differing  wavelength  dependence  which  then 
characterizes  the  solution.  Physically  the  increasing  value  of  S  is  associated  with 
the  specular  lobe,  and  at  a  =  k  the  contribution  (4. 23)  is  present  and  dominates  the 
return  through  its  lower  wavelength  dependence.  Inasmuch  as  these  are  all  first 

+  Vsin(-a)  must  be  interpreted  as  i  Vslna  in  the  formula  for  u(-a). 


94 


THE  UNIVERSITY  OF  MICHIGAN 

3648-2- T 


order  results  only,  they  give  no  Information  about  the  behavior  of  the  field  in  the 
immediate  vicinity  of  a  =  *■,  and  until  such  time  as  more  detailed  results  are  avail¬ 
able,  the  only  method  for  estimating  the  transition  is  to  use  the  physical  optics  solur 
tion  for  a  disc.  We  then  have 


a«_. 


lka 


-2ikacot(x  -  8  )cos  a 


J,  (2ka  sina)  e 


(4.40) 


2  tana  1 

(Crispin  et  al,  1959),  which  can  be  used  for  those  angles  a  near  to  r  at  which  the 
magnitude  of  (4.40)  exceeds  that  of  (4.38). 

In  the  aspect  range  r~9Q  <o  <  x/2  the  Becond  diffracted  wave  doeB  not 
appear,  and  the  geometrical  theory  of  diffraction  then  gives 

S  =  ufcr)  .  (4.41) 

This  is  analogous  to  the  physical  optics  result  (4. 23),  and  each  can  be  written  as 

1  jir/4  -  2ka  cos  (or  -  0J  cosec  0^ 


S  =  A(or,0o)  ^ 


(4.42) 


Xsina 

where  A  is  a  real  quantity,  but  the  values  of  A  are  quite  distinct  both  in  form  and 
in  magnitude.  Unfortunately,  a  study  of  the  experimental  data  does  not  provide  a 
clear-cut  Indication  of  the  formula  which  is  preferable.  Because  of  the  inaccuracy 
associated  with  the  physical  optios  prediction  of  the  return  from  a  shadow  boundary 
(the  base  is  part  of  the  boundary  at  these  aspects)  it  would  seem  that  equation  (4. 41) 
should  be  more  accurate,  and  certainly  it  does  explain  some  of  the  dependence  on 
polarization  observable  in  the  data  (we  shall  later  examine  another  source  of  a 
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polarization-dependent  return  not  predicted  by  either  theory).  On  the  other  hand, 

(4.41)  is  discontinuous*  at  a  =  r-9  (cf  equation  4. 38),  whereas  the  optics  result  is 

o 

continuous,  and  it  is  also  infinite  at  a  *  Q  -  r/2.  This  last  is  due  to  the  specular 
return  from  the  sides  of  the  cone,  and  in  a  neighborhood  of  a  *  r/2  we  have  no 

alternative  but  to  use  (4.29),  replacing  it  by  (4.41)  when  the  magnitudes  of  the  two 
contributions  have  become  comparable. 

In  passing  we  note  that  the  aspect  range  v-0q  <  a  <  »/ 2  is  one  for  which  a 
portion  of  the  conical  surface  is  in  shadow,  and  it  is  now  possible  for  a  creeping 
wave  contribution  to  appear.  This  is  a  polarization-dependent  phenomenon  and  has 
been  discussed  by  Goodrich  et  al  (1958,  1959)  in  connection  with  the  radiation  pat¬ 
tern  of  slots  on  a  cone.  Its  importance  in  the  scattering,  problem  has  not  yet  been 
assessed,  but  it  is  unlikely  to  be  significant. 

Let  us  now  consider  the  cone-sphere.  The  theory  given  by  Keller  (1960)  pre¬ 
dicts  that  for  incidence  within  the  backward  cone  the  scattered  field  is  composed  of 
contributions  from  the  tip  and  from  the  shadow  boundary.  The  first  of  these  can  be 
obtained  using  physical  optics,  and  is  given  in  (4.26).  The  second  is  attributable  to 
creeping  waves  which  have  traversed  the  back  of  the  sphere,  and  since  the  separa¬ 
tion  between  the  cone-sphere  join  and  the  boundary  increases  with  Increasing  fre¬ 
quency,  it  seems  reasonable  to  estimate  this  return  by  reference  to  the  known 

+It  is  also  discontinuous  at  a  *  r/2. 
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sphere  solution.  Following  Keller,  the  appropriate  far  field  amplitude  is  then 

,  ikbf*  +  2oosoooseo(i’-d  )~l  f  -t^kb)1^3  -  ZrJkb)1/3"! 

8  - -±  <kb>4'3  e  L  oJ  0He  H  -Gs.  S 

*“  (4.43)-’ 

where  b  =  a  sec  (t-Q  )  Is  the  radius  of  the  sphere  and  G„,  Gc,  f  T  are  cer- 

tain  parameters  which  are  defined  in  terms  of  Airy  functions  and  their  zeros.  The 

second  term  is  negligible  in  comparison  with  the  first. 

When  the  relative  magnitudes  of  (4.26)  and  (4.43)  are  examined,  it  is  found 

that  for  all  except  the  very  highest  frequencies  the  latter  Is  dominant  In  spite  of  Its 

exponential  decay.  Thus,  for  example,  when  5r/6  and  a  *  0,  equality  between 

the  two  contributions  is  not  reached  until  kb  *  45. 8  and  consequently,  as  kb  in- 

.  2 

creases,  the  back  scattering  cross  section  cr/X  should  display  an  oscillation  of 
increasing  amplitude  which  reaches  a  maximum  when  kb  =  45. 8,  and  thereafter 
decreases.  In  particular,  for  kb  less  than  (about)  15  the  curve  should  be  almost 
monotonic  with  the  cross  section  given  by  the  creeping  waves  alone;  and  for  kb 
greater  than  (about)  80  the  cross  section  should  be  effectively  constant  and  equal  to 
the  tip  contribution. 

These  results  do  not  agree  with  experiment.  Unfortunately  the  most  reliable 
data  is  limited  to  narrow  angle  cones  with  kb  <  14,  and  only  a  few  isola>  jd  measure¬ 
ments  have  been  made  at  higher  frequencies  or  with  cones  of  wider  angle  (0  <  5*76), 

o 

but  on  this  evidence  equations  (4.26)  and  (4.43)  are  Inadequate,  ft  is  clearly  unre- 
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alistic,  of  course,  to  expect  that  the  formulae  will  be  valid  for  kb  of  order  unity 
(even  the  approximation  to  the  creeping  wave  return  is  of  insufficient  accuracy  for 
this  purpose),  but  the  magnitude  and  trends  of  the  measured  values  are  quite  dis¬ 
tinct  from  those  indicated  by  the  theory  even  for  kb  =  0(10).  Throughout  the  fre¬ 
quency  range  5  <  kb  <  14  the  observed  nose-on  cross  section  oscillates  in  a  rela¬ 
tively  uniform  manner  with  a  period  which  is  in  good  agreement  with  that  implied  by 
(4.26)  and  (4.43).  The  average  level  of  the  return,  however,  is  much  greater  than 
predicted,  as  is  the  amplitude  of  oscillation.  This  amplitude  shows  no  significant 
variation  throughout  the  range,  and  certainly  there  is  no  evidence  of  the  expected 
build-up  towards  the  higher  frequencies. 

At  nose-on  incidence  a  considerable  improvement  in  the  theoretical  predic¬ 
tions  can  be  obtained  by  including  a  return  from  the  cone-sphere  join  and  by  using  a 
more  complete  representation  for  the  creeping  wave  contribution.  This  is  essen¬ 
tially  what  has  been  done  by  Kennaugh  and  Moffatt  (1962)  in  the  course  of  applying 
their  "Impulse  approximation"  technique  (Kennaugh  and  Cosgriff,  1958).  From  the 
data  available  in  the  Radiation  Laboratory  it  would  appear  that  the  first  of  these  con¬ 
tributions  can  be  estimated  by  physical  optics  with  an  error  which  is  at  most  20  per¬ 
cent  for  kb  ~  10  .  The  far  field  amplitude  associated  with  the  tip  and  the  Join  is 
then 


_  _  2ikbcos  0  cot  6 

,  2.  2a  oo 

tan  6  -  sec  6  e 
o  o 


(4.44) 
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which  is  believed  to  be  the  correct  high  frequency  solution  for  a  cone-sphere  at 
nose-on  incidence.  For  narrow  angle  cones  the  Inclusion  of  the  second  term  pro¬ 
duces  a  large  increase  in  S,  and  removes  most  of  the  discrepancy  between  theory 
and  experiment  at  the  highest  frequencies  (Id)  ~  80)  for  which  data  is  available.  At 
such  frequencies  the  creeping  wave  contribution  can  be  calculated  with  sufficient 
accuracy  using  (4.43),  but  the  error  increases  with  decreasing  kb.  For  kb  of  ordei 
10  the  error  is  significant,  and  a  more  complete  expression  which  is  correct  to 
within  1  percent  for  kb  >  5  is 

+i»'3_j _ r  i 

W  t  e 


t/3  8g 


'kb f/3  ^ + 2cos a cosec  (r - 0q)} 

n  )  © 


«•« 

(Senior  and  Goodrich,  1963)  where 

0  -  1.01879... 

and  Ai(-0)  =  0.53565... 

At  still  lower  frequencies  even  this  becomes  inadequate  owing  to  the  greater  impor¬ 
tance  of  the  higher  order  waves  which  are  neglected  in  (4. 45),  and  no  simple  formu¬ 
la  for  the  shadow  boundary  return  is  then  available.  On  the  other  hand,  numerical 
values  can  be  obtained  by  subtracting  the  specular  contribution  (Logan,  1960)  from 
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the  complete  far  field  amplitude  computed  from  the  Mie  series.  This  has  been  done 
by  Gent  et  al  (I960)  for  selected  values  of  kb  and.  Implicitly,  by  Kennaugh  and 
Moffatt  (1962). 

hi  spite  of  the  Improvement  In  the  predictions  of  the  nose-on  behavior  result¬ 
ing  from  these  modifications,  there  is  still  a  systematic  discrepancy  between  theory 
and  experiment  which  is  clearly  evident  for  kb  greater  than  (about)  4  and  which 
probably  persists  until  such  time  as  the  creeping  wave  contribution  Is  negligible. 

As  yet  there  Is  no  theoretical  method  for  calculating  the  additional  return,  though 
some  of  Its  properties  can  be  determined  by  detailed  analysis  of  experimental  data. 

For  angles  of  incidence  away  from  nose-on,  the  data  available  for  the  cone- 
sphere  is  less  extensive;  but  on  the  basis  of  those  patterns  which  have  been  taken,  a 
few  general  comments  about  the  back  scattering  cross  section  in  the  resonance  re¬ 
gion  can  be  made.  In  the  Immediate  vicinity  of  a  -  0  the  behavior  depends  on  the 
relative  phases  of  (4.44)  and  (4. 45).  Since  the  polar  diagram  of  the  return  from  the 
Join  is  peaked  at  nose-on  and  falls  to  zero  In  an  angle  comparable  with  the  half  cone 
angle  t-9q,  o/X  has  a  maximum  or  minimum  at  a  *  0  if  (4.44)  and  (4.45)  are  in 
phase  or  out  of  phase  respectively,  changing  over  to  a  minimum  or  maximum 
respectively  at  a  value  of  a  which  is  fractionally  greater  than  r-  .  As  a  In¬ 
creases  beyond  this,  the  cross  section  displays  a  marked  dependence  on  polarization. 

,  2 

With  vertical  polarization  o/X  oscillates  with  a  relatively  small  amplitude 
about  a  level  which  Is  in  good  agreement  with  that  determined  by  (4. 45)  alone,  and 
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this  continues  until  the  side  lobes  of  the  specular  flash  are  reached.  The  latter  can 
be  obtained  from  (4. 29),  and  depending  on  the  phase  relationship  between  (4. 29)  and 
(4.45)  there  may  be  a  few  deep  minima  in  the  aspect  range  for  which  the  two  contri¬ 
butions  are  comparable  in  magnitude.  The  specular  flash  itself  is  in  excellent 
agreement  with  the  scattering  amplitude  shown  in  (4. 30),  but  for  the  side  lobes  be¬ 
yond  a  *  6  -  t/2  it  is  necessary  to  take  into  account  the  presence  of  the  spherical 
cap  in  applying  the  physical  optics  method.  The  results  obtained  when  this  is  done 

are  reasonably  accurate.  For  a  greater  than  (about)  1 J0J  the  cxobb  section  is 

o 

effectively  constant  out  to  180  ,  with  a  value  which  corresponds  quite  closely  to  the 
specular  return  from  the  sphere.  The  appropriate  far  field  amplitude  is  therefore 


-2ikb  cosec  6  (l+sin  9 


{• 


2  sin  6  coso 
o 


1/2 


(4.46) 


(Senior  and  Goodrich,  1963).  Such  discrepancies  as  are  apparent  may  be  due  to  the 
reflection  of  the  creeping  wave  on  the  sphere  at  the  Join  with  the  cone,  but  the  result 
lng  contribution  should  be  small,  and  has  not  been  observed  in  experimental  data 
with  kb  >  8.  We  note  in  passing  that  Moffett  (1962)  has  proposed  a  modification  of 
the  physical  optics  Integral  for  use  in  predicting  the  oross  section  at  or  *  180°.  The 
agreement  with  experiment  is  probably  within  the  experimental  error.  On  the  other 
hand,  the  data  obtained  by  the  Radiation  Laboratory  does  not  reveal  the  need  for  any 
modification  (or  addition)  to  (4.46). 
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With  horizontal  polarization  the  noae-on  peak  or  trough  is  somewhat  nar- 

/  2 

rower  than  above,  and  beyond  the  turning  point  at  oA  &lls  away  rapid¬ 

ly.  Apart  from  a  few  large  but  isolated  peaks  the  cross  section  out  to  the  side  lobes 
of  the  specular  flash  is  many  db  below  that  for  vertical  polarization,  and  this  is 
explainable  by  the  absence  of  the  dominant  creeping  wave  return  (the  cone  sides 
prevent  any  such  return  from  reaching  the  receiver  without  first  undergoing  a  scat¬ 
tering  at  the  tip).  Because  of  the  striking  difference  in  behavior  at  the  two  polari¬ 
zations  it  would  not  appear  that  tbs  cylindrical  current  method  (which  is,  in  prac¬ 
tice,  independent  of  polarization)  is  applicable  to  this  aspect  range,  though  its  use 
has  been  advocated  by  Moffatt  (1962). 

Through  the  region  of  the  specular  flash  and  its  side  lobes,  the  cross  section 
shows  little  or  no  variation  with  polarization,  and  once  again  the  physical  optics 
method  can  be  employed  in  its  prediction,  but  beyond  a  ~  100°  (approx. )  an  increa¬ 
sing  oscillation  sets  in,  and  out  to  180°  the  behavior  is  quite  distinct  from  that  ob¬ 
served  with  vertical  polarization.  The  average  level  is  somewhat  greater  than  that 
indicated  by  the  sphere  return  alone,  and  the  oscillation  increases  in  amplitude 
reaching  a  maximum  some  few  degrees  short  of  rear-on.  At  or  *  180°  the  cross 
section  varies  from  local  maxima  to  local  minima  as  a  function  of  frequency  and  its 
value  here  is  Identical  to  that  with  vertical  polarization. 

It  is  believed  that  these  oscillations  are  attributable  to  travelling  waves  ex¬ 
cited  on  the  surface  of  the  body.  Such  waves  can  produce  a  significant  contribution 
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to  the  back  scattering  cross  section  only  with  horizontal  polarization,  and  the  theory 
associated  with  them  predicts  an  aspect- sensitive  return  which  is  similar  to  that 
found  at  near  rear-on  aspects  with  the  cone-sphere,  and  in  the  aspect  range  between 
a  -  *-©o  and  the  side  lobes  of  the  specular  flash  for  a  cone  with  almost  any  termi¬ 
nation. 

When  a  long  thin  body  is  illuminated  at  near  not  .-on  incidence  the  initial 
scattering  is  predominantly  in  the  forward  direction.  A  portion  of  the  energy  is 
guided  along  the  surface,  leading  to  a  current  distribution  which  is  identifiable  with 
a  wave  travelling  down  the  body.  When  this  reaches  the  further  end  a  fraction  of  it 
is  reflected  back  in  the  direction  from  which  it  came,  and  ultimately  radiates  in  the 
back  scattering  direction.  As  a  result  the  observed  signal  appears  to  originate  at 
(or  near)  the  rear  of  the  body,  and  suggested  to  Peters  (1958)  and  Belkina  (1957)  that 
this  type  of  structure  should  act  as  a  travelling  wave  antenna.  By  likening  it  to  a 
long  thin  wire,  Peters  obtained  an  expression  for  the  far  field  amplitude  in  the  form 

S  =  yF(O)  (4.47) 

where  7  is  the  (complex)  current  reflection  coefficient  and  F(0)  is  a  pattern  factor 
which  is  a  function  of  the  length  of  the  wire  in  wavelengths  and  of  the  effective  phase 
velocity  along  the  surface. 

It  is  obvious  that  the  method  is  only  as  accurate  as  is  the  simulation  by  the 
wire,  and  is  certainly  limited  in  application  to  bodies  whose  length  is  large  in  com¬ 
parison  with  both  the  wavelength  and  the  width.  Because  of  the  basic  model  the 
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results  are  markedly  dependent  on  the  polarization,  and  the  travelling  waves  provide 
a  significant  contribution  to  the  oross  section  only  at  angles  away  from  nose-on 
(F  s  0  for  6  ~  0)  with  horizontal  polarization.  With  any  actual  body  having  signi¬ 
ficant  width,  a  small  contribution  with  the  other  polarization  Is  to  be  expected,  but 
cannot  be  estimated  using  Peters'  theory. 

In  spite  of  the  inaccuracies  inherent  in  the  model,  the  method  has  proved  ex¬ 
tremely  useful  in  practice.  Applications  to  the  finite  cone  have  been  given  by  Cris¬ 
pin  et  al  (1959)  and  Dawson  (1960b),  and  to  the  cone-sphere  by  Gent  et  al  (1960),  and 
for  such  purposes  it  is  generally  sufficient  to  assume  that  the  phase  velocity  along 
the  surface  is  the  same  as  that  in  free  space.  The  expression  for  F(0)  then  reduces 
to 

cot2  |  2  2  0 

F(9>  ‘  Iml-  0.4228  (k/*U‘  2>  (4  48> 

where  /  is  the  overall  length  of  the  body  (or  the  distance  from  the  nose  to  the 
reflection  point,  if  this  is  different) .  The  maximum  in  the  pattern  occurs  at 

0  ~  49.35ft.//)1/2  degrees,  (4.49) 

which  decreases  as  kjL  increases.  In  a  neighborhood  of  this  angle  the  pattern  is 
2 

similar  to  (sin x/x)  in  shape.  The  maximum  value  of  S  is 


0.725  ykl 

max*  ~jLzkl  -  0.4228 


(4.50) 
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and  the  Important  role  played  by  the  reflection  coefficient  la  now  apparent.  It  la 

therefore  unfortunate  that  no  method  is  yet  available  for  computing  y  In  terms  of  the 

shape  parameters  of  the  body.  By  analysing  a  large  volume  of  experimental  data, 

2 

Dawson  (1960a)  has  suggested  that  y  is  proportional  to  sin  6,  with  a  constant  of 
proportionality  which  decreases  from  unity  if  the  termination  is  sharp,  to  zero  when 
the  radii  of  curvature  at  the  termination  are  several  wavelengths  or  greater.  There 
is  no  theoretical  support  for  this  suggestion  and  the  dependence  on  6  appears  ques¬ 
tionable.  Other  (quite  different)  estimates  of  y  have  been  given  by  Siegel  (1959), 
but  even  without  a  precise  knowledge  certain  general  statements  can  be  made. 

n  the  termination  is  sharp,  y  is  almost  certainly  independent  of  X  to  the 
first  order:  then  decreases  with  decreasing  X  at  an  extremely  slow  rate,  and 

in  most  practical  cases  (4. 50)  will  be  comparable  in  magnitude  to  the  nose-on  return 
from  a  flat  base  (equation  4.32).  It  will  therefore  be  the  dominant  feature  of  the  re¬ 
turn  in  the  aspect  range  between  the  nose-on  and  specular  regions  when  the  polariza¬ 
tion  is  horizontal,  and  the  large  peaks  which  are  observed  are  in  approximately  the 
positions  predicted  by  the  theory.  It  will  also  be  significant  at  rear-on  aspects  for 
a  cone-sphere,  since  the  reflection  point  is  now  the  apex. 

If  the  rear  of  the  body  is  smooth  the  reflection  coefficient  will  be  much  smal¬ 
ler,  and  depending  on  the  precise  shape  of  the  body  the  reflection  may  occur  at  some 
point  where  a  higher  derivative  of  the  surface  profile  has  a  sharp  discontinuity. 

This  is  believed  to  be  the  case  with  the  cone-sphere,  and  the  experimental  data  on 
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the  scattering  cross  section  at  near  nose-on  aspects  is  suggestive  of  a  reflection  at 
the  cone-sphere  join,  with  y  proportional  to  X  (or  possibly  some  higher  power). 
Since  the  nose-on  return  also  has  this  wavelength  dependence,  the  relative  impor¬ 
tance  of  the  travelling  waves  is  unchanged . 
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EXPERIMENTAL  DATA 

The  Interest  in  cone- like  structures  has  Increased  considerably  during  the 
last  few  years  and  a  wealth  of  experimental  data  is  now  extant.  Unfortunately  a  good 
deal  of  this  still  exists  only  in  the  files  of  the  laboratories  responsible  and  will 
doubtless  remain  so  until  some  central  organization  can  undertake  (or  delegate)  the 
compilation  task.  Because  of  the  needs  of  a  particular  laboratory  at  a  given  time, 
the  data  which  they  obtained  may  have  been  restricted  to  a  few  isolated  cases  and  is 
in  that  sense  "unpublishable'',  but  when  viewed  in  conjunction  with  comparable  data 
obtained  elsewhere  it  could  be  extremely  valuable  in  filling  in  the  "gaps"  which  exist. 

The  following  listing  is  the  result  of  a  search  through  the  open  and  classified 
literature,  together  with  such  unpublished  data  as  the  authors  are  personally  aware 
of.  It  is  certain  that  the  list  is  not  complete  and  it  may  even  contain  only  a  small 
fraction  of  the  measurements  which  have  been  made.  We  would  be  grateful  for  infor¬ 
mation  about  any  data  which  we  have  unintentionally  omitted. 

In  almost  every  case  the  quantity  which  has  been  measured  is  the  back  scat¬ 
tered  cross  section,  and  attention  will  be  confined  to  this;  little  is  known  about  the 
phase  or  about  the  bistatlc  behavior.  The  flat  backed  cone  is  considered  first,  fol¬ 
lowed  by  the  cone-sphere  (referred  to  as  the  "carrot”  in  some  British  publications) 
and,  in  the  third  section,  by  a  brief  resume  of  the  data  for  cones  with  other  types  of 
terminations. 
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5. 1  Flat- Backed  Cone 

Let  a  denote  the  radius  of  the  base  and  t-0q  the  half  cone  angle. 

Using  cw  equipment  operating  at  3.02,  3.8,  9.3  and  9.9  KMo,  Brysk  et  al 
(1959)  have  measured  the  nose-on  cross  section  in  the  following  oases: 

ka  =  0.118  to  3.0  (25  values) 
ka  =  0.143  to  3.0  (21  values) 
ka  =  0. 165  to  3. 0  (26  values) 
ka  =  0.327  to  3.2  (20  values) 
ka  =  0.490  to  4.8  (20  values) 


t-8  =  7.5" 
o 

x-8  m  9.6° 
o 

t-8  *  12° 
o 

r-8  =  45° 
o 


t-8  =  60 
o 

Data  was  also  obtained  for 


t-0  *  4°  :  ka  =  0.187  to  1.05  (14  values) 

o 

using  the  same  equipment  at  8. 465  and  9. 9  KMc.  The  presentation  is  in  the  form 
.  2 

of  graphs,  o/ra  versus  ka,  on  a  log- log  scale. 

Keys  and  Primlch  (1959a)  have  used  two  cw  systems  operating  at  8.75  and 

35  KMc  to  obtain  the  nose-on  cross  section  for  cones  of  half  angle  t-8  =  4°,  20° 

o 

and  60°  with  ka  varying  from  0. 5  to  9. 5  (approx. ).  The  data  is  presented  as 
2 

graphs  of  o/X  versus  2a/X  on  a  log-log  scale.  The  results  for  the  widest  angle 
cone  are  partially  reproduced  in  Keys  and  Primlch  (1959b),  and  corresponding  values 
for  t-0q  =  7. 5°  and  30°  are  also  provided. 

Probably  the  most  complete  set  of  data  is  that  given  in  Keys  and  Primlch 
(1959c  and  d).  A  pulsed  system  was  used  operating  at  35  KMc  with  a  pulse  length  of 
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10  millimicroseconds.  The  back  scattering  cross  section  cr/X  is  shown  as  a  func¬ 
tion  of  aspect,  0  to  360°,  for  vertical  and  horizontal  polarization  and  for  the  follow¬ 
ing  values  of  ka: 

ka  -  0.933,  1.40,  1.86,  2.33,  3.08,  3.80, 

4.56,  5.31,  6.06,  6.79,  7.45,  8.29, 

In  the  first  of  the  above  reports  the  cone  angles  treated  have 

r-Q  *  4°,  7.5°,  9.6°,  12°,  15°  and  20° 
o 

and  in  the  second 

w-eQ  -  30°,  37.5°,  45°,  52.5°  and  60°. 

The  overall  accuracy  is  stated  as  1 2  db. 

If  the  nose-on  cross  sections  are  read  off  from  these  patterns  and  converted 
2 

to  values  of  o,  t a  ,  the  points  shown  in  Figs.  10  through  20  are  obtained.  The 
curves  have  been  calculated  using  equations  (4. 7)  and  (4.37)  (Siegel's  modified  Ray¬ 
leigh  theory  and  Keller's  second  order  theory),  and  are  in  good  agreement  with  the 
data.  Four  typical  aspect  patterns  are  reproduced  in  Figs.  21  through  24  and  in  the 
following  four  figures  the  measured  cross  sections  for  the  specular  flash  are  plotted*] 
as  functions  of  ka  for  jr-0Q  =  4°,  15°,  30°  and  60°.  The  straight  lines  are  the 
physical  optics  predictions  (see  equation  4. 30).  In  most  cases  the  discrepancy  is  no 
more  than  a  factor  2  at  most  and  because  of  the  nature  of  the  theoretical  method  the 

+The  values  for  the  left  and  right  flashes  are  both  presented  if  they  differ  by  more 
than  1  db. 


4.18 

9.02. 
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agreement  improves  with  increasing  t-0q  or  ka.  We  note  in  passing  that  the  cylin¬ 
drical  current  prediction  (also  plotted  in  Figs.  25-28)  is  smaller  than  the  physical 
4 

optics  by  a  factor  see  Qq  in  the  cross  section.  The  factor  increases  from  1 . 0096 
o  o 

when  r-Q  =4  to  4  when  r-Q  *  60  and  on  this  basis  it  would  appear  that  the 
o  o 

"approximated  form"  of  the  cylindrical  current  method  is  somewhat  inferior  to 
physical  optics. 

For  a  cone  of  half  angle  r-Q  =  7. 5°  a  series  of  measurements  of  the  nose- 

o 

on  and  base-on  cross  sections  have  been  made  by  Olte  and  Silver  (1959)  using  a  cw 
ground  plane  system  at  9.33  and  9. 331  KMc.  The  base  radius  was  changed  by 
adding  "discs"  to  one  of  several  different  basic  cones,  giving  a  coverage  of  the  range 

ka  *  1.90  to  6.20 

with  a  total  of  15  values  shown  for  nose-on  and  14  for  base-on.  The  results  are  pre¬ 
sented  in  the  form  of  graphs  of  o/cq  (in  db)  versus  L/\  on  a  linear  scale,  where  L 
is  the  length  of  the  cone  and  oq  is  the  scattering  cross  section  of  a  6  in.  diameter 
reference  sphere. 

The  same  equipment  was  earlier  employed  by  Shostak  and  Angelakos  (1957)  to 
measure  the  scattering  pattern  as  a  function  of  aspect  for  a  cone  of  half-angle  7. 5°. 
At  that  time,  however,  the  accuracy  of  the  system  was  poor  and  the  experiments 
were  later  repeated  by  August  and  Angelakos  (1960).  For  w -9  *  7. 5°  and 
ka  *  1. 53,  3. 14  and  6. 18  the  measured  nose-on  and  base-on  cross  sections  are 
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listed  and  in  each  case  the  aspeot  pattern,  0  to  180°,  for  vertical  polarization  is 
shown.  Measurements  on  these  same  models  were  also  reported  by  Honda  et  al 
(1959). 

Additional  results  for  cones  with  «=  7. 5°  have  been  obtained  by  the  Cor¬ 
nell  Aeronautical  Laboratory  and  Microwave  Radiation  Co. ,  Inc. ,  and  the  measured 
values  of  the  nose-on  cross  sections  were  communicated  to  K.  M.  Siegel.  No  infor¬ 
mation  about  the  equipment  is  available,  but  for  ease  of  reference  the  data  is  listed 
below: 


o/X*  (in  db) 

ka 

Hbr.  Pol. 

Vert.  Pol. 

Organization 

38.7 

19.0 

19.0 

Cornell  Aeronautical  Lab. 

23.7 

13.4 

13.4 

7.4 

7.5 

7.2 

26.4 

21.9 

17.3 

Microwave  Rad.  Co. ,  Inc. 

6.6 

1.2 

2.1 

3.3 

-6.0 

-5.4 

1.64 

-2.8 

-3.1 

The  only  remaining  measurement  of  which  we  are  aware  was  carried  out  at 

the  Royal  Aircraft  Establishment  in  England  using  a  pulsed  system  at  X  band 

(X  -  3.2  om).  For  r-9  ■  15°  and  ka  *  9.016  the  scattering  pattern  as  a  function 

o 

of  aspect,  0  to  120°  has  been  published  by  Dawson  and  Turner  (1960).  The  results 
for  both  polarizations  are  given,  and  in  Dawson  (1960b)  the  aspeot  range  is  extended 
out  to  180°. 
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5.2  The  Cone-Sphere 

Much  of  the  original  data  on  the  cone-sphere  has  appeared  only  In  classified 

publications  and  this  Is  true,  for  example,  of  the  work  carried  out  at  the  Royal 

Radar  Establishment  using  a  cw  doppler  system  at  9. 37  KMc  (X  *  3. 20  om).  Gent 

et  al  (1960)  have  reported  measurements  on  50  similar  cone-sphere  models  with 

j-8  *  15°  and 
o 

b/X  *  1.00  (0.01)  1.06,  1.08,  1.105,  1.125(0.025)2.125 
where  b  Is  the  sphere  radius.  With  each  model  the  back  scattered  cross  section 
was  recorded  for  horizontal  polarization  at  1°  Increments  in  the  aspect  angle  In  the 
range  t  5°  about  the  nominal  nose-on  position  and  from  the  symmetry  of  the  result¬ 
ing  curves  the  precise  value  of  the  nose-on  cross  section  was  determined.  In  six 
cases  (b/X  *  1.01,  1.08,  1.125,  1.925,  2.00  and  2.075)  the  patterns  for 
-2°  <  a  <  20°  (or  greater)  are  presented,  together  with  the  complete  aspect  pattern, 
0  to  360°,  for  the  model  having  b/X  *  1 . 02 . 

During  tim  past  year,  however,  data  on  the  cone-sphere  has  started  appear¬ 
ing  In  the  open  literature.  The  first  such  release  was  by  Kennaugh  and  Moffatt  (1962] 
who  used  a  cw  X-band  radar  to  measure  the  nose-on  cross  section  of  a  series  of 
cone-spheres  having 

w-9  *  15°  and  kb  -  1  to  7  (39  values), 
o 

The  spacing  of  the  points  provides  a  reasonably  complete  coverage  of  this  range  In 
kb  and  serves  to  extend  the  data  of  Gent  et  al  (1960)  almost  to  the  edge  of  the 
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i  2 

Rayleigh  region.  The  results  are  presented  in  the  form  of  a  graph  of  o/rb  in  db 

versus  kb  on  a  linear  scale.  A  larger  version  of  the  graph  has  been  provided  by 

Moffett  (1962).  These  points  have  been  renormalized  and  appear  In  Fig.  29. 

The  most  extensive  data  so  far  obtained  has  recently  been  reported  by  Blore 

and  Royer  (1962).  Using  the  8. 6  mm  pulse  equipment  of  the  Defence  Research  Boarc 

in  Canada  the  following  nose-on  cross  sections  have  been  measured: 

x-9  *  7.5°  kb  =  0.4  to  2.8  (30 values,  approx.) 
o 

t-Q  -  15°  kb  *  0.2  to  9.6  (105 values,  approx.) 

o 

»-0q  *  20°  kb  *  0.5  to  9.5  (80  values,  approx.) 

x-0q  *  30°  kb  *  0.2  to  9.6  (105  values,  approx.) 

x-0  ■  37. 5°  kb  *  0.2  to  5. 7  (85 values,  approx.) 

o 

All  values  are  the  average  of  at  least  two  measurements  with  both  horizontal  and 

vertical  polarization,  and  are  believed  accurate  to  t  0. 5  db.  The  results  are  pre- 

sented  In  the  form  of  graphs  of  c/X  in  db  versus  26/X  on  a  linear  scale,  and  In 

Figs.  29  and  30  we  have  reproduced  the  curves  for  two  of  the  cone-spheres. 

Isolated  measurements  on  single  models  have  also  been  carried  out  at  the 

Royal  Aircraft  Establishment  In  England  using  an  X-band  pulse  equipment,  and  the 

o 

data  Is  available  In  their  reports.  For  the  cone-sphere  having  t-0q  *  15  and 
kb  ■  44. 9  Dawson  and  Turner  (1960)  give  the  aspect  patterns  (0°  to  80°)  for  horizon¬ 
tal  polarization,  and  the  analogous  results  for  kb  =  6.28  have  been  presented  by 
Dawson  (1960b).  With  this  same  cone  angle,  but  kb  *  18. 70,  43.89  and  89.77, 
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Turner  (1960)  has  provided  aspect  patterns  out  to  70°  for  vertical  polarisation,  and 
also  shows  the  results  for  a  cone-sphere  with  r-0Q  s  25°  and  ka  *  30. 0. 

The  above  is  a  complete  listing  of  the  published  data  on  the  cone-sphere 
known  to  the  present  authors,  but  it  Is  certain  that  many  of  the  scattering  measure¬ 
ments  for  this  shape  have  not  yet  been  reported  In  the  literature.  As  an  Illustration 

o 

of  this  fact,  the  following  measurements  on  a  cone-sphere  having  t-0q  *  12. 5  and 

b*4. 5187  cm  have  been  carried  out  by  the  Radiation  Laboratory,  and  this  is  only 

one  of  the  many  cone-sphere  models  whloh  have  been  Investigated: 

Aspect  patterns  0°  to  360°: 

2 . 67  (0. 10)  3.07  KMc,  horizontal  polarization 

8. 58  and  9. 375  KMc,  horizontal  and  vertical  polarization 

8.5,  9.3,  10.0  and  31 . 97  KMo,  vertical  polarization 

Aspect  patterns  0°  to  50°: 

7.98,  8.1,  8.34,  8.43,  8.58,  8.8,  8.93,  9.0, 

9.15,  9.3,  9.6,  9.8,  9.954,  10.11,  10.43,  10.97, 
and  11.97  KMc,  vertlo&l  polarization. 

At  each  of  these  frequencies  16  or  more  Independent  measurements  of  the  noee-on 
cross  section  have  been  made  to  ensure  the  most  accurate  values  possible  over  a 
complete  cycle  of  the  osolllation  described  In  Section  4. 3.  All  measurements  were 
made  with  a  cw  system.  It  Is  hoped  to  publish  at  least  some  of  this  data  in  the  near 
future. 
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5.3  Other  Ter™<na*<™m 

For  terminations  other  than  a  flat  back  or  a  sphere  with  no  discontinuity  in 
the  first  derivative  at  the  Join,  a  comparison  of  the  experimental  data  is  more  diffi¬ 
cult  owing  to  the  variations  in  the  structural  features  of  the  individual  models. 
Nevertheless,  it  may  be  helpful  to  mention  some  of  the  data  which  is  available  for 
more  general  oonical  shapes. 

For  two  flat-baoked  cooes  mounted  base  to  base  with  w-d  *  20°,  Blore  and 

o 

Royer  (1962)  have  measured  the  nose-on  cross  section  for 

ka  *  0.3  to  9.3  (60  values,  approx.) 

where  a  is  the  radius  at  the  join.  These  same  authors  have  also  given  data  for 

doubly-rounded  cones  (the  above  shape,  but  with  the  Junction  of  the  two  cones  roundec 

off)  with  w-6  -  7. 5°,  20°  and  37. 5°. 
o 

Four  different  rounded  flat-booked  cones  have  been  investigated  by  Turner 
(1960)  and  for  each  the  X-band  aspect  pattern,  0°  to  70°,  for  vertioal  polarization  is 
presented.  If  a  denotes  the  maximum  radius  and  6  is  the  radius  of  curvature  at 
the  Join  of  the  cone  and  the  base,  the  models  studied  have  the  following  specifications 


w-6  *  15° 

0 

ka  *  43.89 

6/X  * 

3.17 

w-6  -  15° 

0 

ka  *  44.89 

<yx. 

3.17 

w-6  -  15° 

0 

ka  -  89.77 

S/\  * 

3.17 

w-6  -  15° 

0 

ka  -  89.77 

S/\  - 

6.35 
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FIGURE  10:  NOSE-ON  RADAR  CROSS  SECTION  OF  FLAT  BACKED  CONE, 


Modified  Rayleigh  Theory  (Siegel  1959} 
Keller's  Second  Order  Theory 
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FIGURE  12:  NOSE-ON  RADAR  CROSS  SECTION  OF  FLAT  BACKED  CONE,  jt-0  =9.6 


Modified  Rayleigh  Theory  (Siegel  1959) 
Keller's  Second  Order  Theory 


FIGURE  13:  NOSE-ON  RADAR  CROSS  SECTION  OF  FLAT  BACKED  CONE, 


FIGURE  14:  NOSE-ON  RADAR  CROSS  SECTION  OF  FLAT  BACKED  CONE, 


Modified  Rayleigh  Theory  (Siegel  1959) 
Keller's  Second  Order  Theory 
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FIGURE  15:  NOSE-OX  RADAR  CROSS  SECTION  OF  FLAT  BACKED  CONE, 
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FIGURE  16:  NOSE-ON  RADAR  CROSS  SECTION  OF  FLAT  BACKED  CONE,  w-0  =  30 
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FIGURE  18:  NOSE-ON  RADAR  CROSS  SECTION  OF  FLAT  BACKED  CONE 
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FIGURE  19:  NOSE-ON  RADAR  CROSS  SECTION  OF  FLAT  BACKED  CONE,  *-0  =52.5 
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FIGURE  20:  NOSE-ON  RADAR  CROSS  SECTION  OF  FLAT  BACKED  CONE 
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FIGURE  21:  BACK  SCATTERED  RADAR  CROSS  SECTION  OF  FLAT  BACKED  CONE 

*  -0  =15°,  ka  *  3.08,  Horizontal  Polarization  (Keys  and  Primich,  1959c) 
o 


FIGURE  22:  BACK  SCATTERED  RADAR  CROSS  SECTION  OF  FLAT  BACKED  CONE 

j r-0  =15°,  ka  =  3.08.  Vertical  Polarization  (Keys  and  Primich,  1959c) 
o 
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[GURE  23:  BACK  SCATTERED  RADAR  CROSS  SECTION  OF  FLAT  BACKED  CONE 
T  -  dQ  =  15°,  ka  =  9. 02,  Horizontal  Polarization  (Keys  and  Prlmich,  1959c) 


FIGURE  24:  BACK  SCATTERED  RADAR  CROSS  SECTION  OF  FLAT  BACKED  CONE 
jr  -  0  =15°,  ka  =  9. 02,  Vertical  Polarization  (Keys  and  Primich,  1959c) 
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FIGURE  25:  RADAR  CROSS  SECTION  OF  FLAT  BACKED  FINITE  CONE 
SPECULAR  FLASH  o(8  -t/2),  t-B  -4° 
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FIGURE  26:  RADAR  CROSS  SECTION  OF  FLAT  BACKED  FINITE  CONE 
SPECULAR  FLASH  o{9 -v/2),  t-6  -7.5° 
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FIGURE  27:  RADAR  CROSS  SECTION  OF  FLAT  BACKED  FINITE  CONE  - 
SPECULAR  FLASH  <r(9-w/2),  t-9  -15° 
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FIGURE  28:  RADAR  CROSS  SECTION  OF  FLAT  BACKED  FINITE  CONE  - 
SPECULAR  FLASH  o(0-r/2),  T-0  -30° 
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FIGURE  29 
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